Overview of the spin structure function g\ at arbitrary x and Q 2 



B.I. Ermolaev 

Ioffe Physico-Technical Institute, 194021 St. Petersburg, Russia 
M. Greco 

Department of Physics and INFN, University Rome III, Rome, Italy 

S.I. Troyan 

St. Petersburg Institute of Nuclear Physics, 188300 Gatchina, Russia 

In the present paper we summarize our results on the structure function gi and present explicit 
expressions for the non-singlet and singlet components of g\ which can be used at arbitrary x 
and Q 2 . These expressions combine the well-known DGLAP-results for the anomalous dimensions 
and coefficient functions with the total resummation of the leading logarithmic contributions and 
the shift of Q 2 — » Q 2 + fi 2 , with (i/Aqcd ~ 10 (~ 55) for the non-singlet (singlet) components 
, of gi respectively. In contrast to DGLAP, these expressions do not require the introduction of 

singular parameterizations for the initial parton densities. We also apply our results to describe the 
£ — ' experimental data in the kinematic regions beyond the reach of DGLAP. 

PACS numbers: 12.38.Cy 

II ■ I. INTRODUCTION 

As it is well-known, the spin structure function g\ is introduced through the following conventional parametrization 
of the spin-dependent part Wf^ n of the hadronic tensor of the Deep Inelastic lepton- hadron Scattering: 
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qq \ where we have used the standard notations: P is the hadron momentum, Mh is the hadron mass, S is the hadron spin, 
^vq ■ q is the virtual photon momentum. Traditionally, — q 2 = Q 2 > and x = Q 2 /2Pq. The scalar functions gi{x, Q 2 ) and 
g2{x,Q 2 ) are called the spin structure functions. Both of them contribute to the asymmetry between the DIS cross 
sections when the lepton and hadron spins are antiparallcl and parallel. In particular, gi describes such asymmetry 
when both spins are longitudinal, i.e. they lie in the plane formed by P and q. Obviously, in order to calculate the 
structure functions gi(x, Q 2 ) and g2(x,Q 2 ) in the framework of QCD, one should know the QCD behaviour at large 
and small momenta of virtual particles, i.e. one should be able to account for the perturbative and non-perturbative 
effects. At present this is impossible, so the standard description involves the factorization hypothesis: W spm [iv is 
represented as a convolution: 



Wp n = W% ®-* q + W% ® * 9 (2) 

where q and ty g are the probabilities to find a polarized quark or gluon in the polarized hadron, while describe 
the DIS of the quarks and gluons. There is no rigorous proof of this factorization, especially at small x, in the 
literature. However, discussing this point is beyond the scope of our paper. Also, there is no model-independent 
theoretical description of the probabilities "Jg and ^ g in the literature because QCD at small momenta is not known. 
On the contrary, W^ v and W^ v can be calculated with the methods of Pertubative QCD, by summing the contributions 
of the involved Feynman graphs. So, the standard procedure is to replace $ ' q and \V g by the initial parton densities 
Sq and Sg. Both of them are found by fitting the experimental data at large x and not very large momenta Q 2 
(Q 2 = [i 2 - 1 GeV 2 ). Therefore, 

.91 = g\ ®Sq + g\ ® Sg. (3) 
It is well-known that in the Born approximation 

gf™ - (e 2 g /2)S(l -x)®5q (4) 
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where e q is the electric charge of the quark interacting with the virtual photon. Accounting for the QCD radiative 
corrections to gf ° rn and other DIS structure functions, especially by trying to perform a complete resummation of 
the corrections, has been the subject of great interest in recent years. Surely, such resummation cannot be performed 
precisely, so it would be important to resum, in the first place, the most essential corrections. They arc different for 
different values of x and Q 2 . For example for describing g\ in the region of x ~ 1 and large Q 2 , the contributions 
~ In x are negligibly small compared to ln fc (Q 2 //i 2 ). In contrast, ln fe x becomes quite important at x -C 1 and should 
be accounted for. 

The goal of obtaining an universal description of the structure function gi, which could be used for arbitrary x 
and arbitrary Q 2 would be appealing both for theorists and experimentalists. Most generally, one encounters various 
kinematic regions where the DIS structure functions have been thoroughly studied. The first kinematic region is the 
so-called hard region A of large x and large Q 2 : 

A: w > Q 2 > ^ 2 , x < 1 (5) 

where w = 2pq and /i 2 is the starting point of the Q 2 -evolution. Usually, the value of /x is chosen w 1 GeV or 
so. Through the paper we use the standard notations: q is the virtual photon momentum, p is the initial parton 
momentum, and x = Q 2 /w. The region A was described first by the LO DGLAP evolution equations obtained in 
Ref. [J: 



dAq 
dt 

dAg 
dt 



P qq ®Aq + P qg ® Ag , (6) 
P gq ®Aq + P gg <g> Ag 



where we have used the standard notation t = ln(Q 2 //i 2 ) and Pik (with i,k — q,g) are the splitting functions. Aq and 
Ag are the evolved (with respect to Q 2 ) parton distributions. The splitting functions Pik in the DGLAP evolution 
equations ^ include the QCD coupling a s . In order to account for the running a s -effects, one should define the 
argument of a s . The DGLAP -prescription is 

a s =a s {Q 2 ). (7) 

Through the paper we will address the parametrization of a s in Eq. ([7|) as the DGLAP -parametrization. The DGLAP 
Eqs. (J6j) describe the Q 2 -evolution of the parton distributions from Q 2 — [i 2 , with n ~ 1 GeV, to larger Q 2 . When 
general solutions to Eqs. © are obtained, one needs to specify appropriate initial conditions. Conventionally, the 
initial conditions to Eqs. ([6]) are 

Aq\ t=0 = 5q , Ag\ t=0 = Sg , (8) 

with 6q, 5g being called the initial parton densities. They are found by fitting the experimental data . After Aq and 
Ag have been fixed, the DIS structure functions, including gi, are found by convoluting them with the coefficient 
functions C q , C g : 

gi(x, Q 2 ) = C q (x/y) ® Aq(y, Q 2 ) + C g (x/y) ® Ag(y, Q 2 ) . (9) 

The Mellin transformations of Pik are called the anomalous dimensions. The splitting functions P^ and coefficient 
functions Ck for the unpolarized DIS were calculated with LO accuracy in Ref. The LO expressions for Pik 
and Ck for the polarized DIS were obtained in Ref. ,2|. Later, the LO expressions of Refs [l], 0] for Pik and Ck 
were complemented by the NLO results 0|. A detailed review on that subject can be found in Ref. [![. From pure 
theoretical grounds, this approach should not be used outside the region A. However, introducing special fits]!, [|[ for 
the initial parton densities DGLAP has been extended to the region B of large Q 2 and small x: 

B: w > Q 2 > n 2 , x < 1 . (10) 

Indeed the parameterizations for Sq, Sg of Ref. contain singular factors x~ a , and used in Eq. ((HJ, they provide 

gi with a fast growth at small x. As a result, combining the LO evolution equations of Ref. [l( and NLO DGLAP 
results of Ref. [3[ with the standard fits of Ref. [H, Q it has been possible to describe the available experimental data 
on gi in regions A and B, i.e. for large Q 2 and arbitrary x. In the present paper we refer to this as the Standard 
Approach (SA). In addition to regions A and B, there are two more interesting kinematic regions: 

C: 0<Q 2 </i 2 , x«l, (11) 
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D: < Q 2 < fi 2 , x < 1 . (12) 

Besides a purely theoretical interest, the knowledge of g\ in the regions C and D is needed because they correspond 
to the kinematic region investigated experimentally by the COMPASS collaboration. Obviously, the regions C and 
D are beyond the reach of SA. Strictly speaking, the same could be said about the region B: In fact the expressions 
for Pik are obtained (see Refs. P, 0| for detail) under the assumption of the ordering 

H 2 < k 2 ± < k 2 ± < ... < Q 2 (13) 

where ki ± are the transverse momenta of virtual ladder partons and they are numbered from the bottom of the 
ladders to the top. Once this ordering is kept one is led inevitably to neglect the double-logarithmic (DL) contributions 
~ a s ln 2 (l/a;) and other contributions independent of Q 2 . Such contributions are small in the region A where they 
are correctly neglected in the SA. However, they become essential in the region B. In order to account for them, the 
DGLAP-ordering of Eq. (fT"3|) should be replaced by the other ordering: 

m2 < fcjLk < < ... < w ( i4) 

Pl P2 

where /3j are the longitudinal Sudakov variables 1 for the virtual parton momenta kj as follows: 

kj = -ctj(q + xp) + [3jp + kjx ■ (15) 

In order to account for such logarithmic contributions, Eq. (|14| should be implemented by the ordering for /3j: 

1 > 0i > fo > ■ ■ > fi 2 /w . (16) 

This ordering does not exist in DGLAP because in this approach /3j > x ~ 1. The ordering (|14I16[) was first introduced 
in Ref. Q in the context of QED but it applies in QCD as well. Replacing the ordering lfl3|l by Eqs. (fT4l [T6|) makes 
possible to sum up all DL contributions, regardless of their argument, to all orders in a S7 i.e. to perform calculations 
in the double- logarithmic approximation (DLA). Explicit expressions for g\ in DLA were obtained in Ref. [1 01 ] . The 
drawback of those expressions is that a s is kept fixed at an unknown scale. The effect of running a s were taken 
into account in Ref. [1JJ. The parametrization of a s in Refs. [Till di ffers from the DGLAP- parametrization. The 
theoretical grounds for this new parametrization were given in Ref. [12J and a numerical comparison with the standard 
parameterizations can be found in Ref. [HI . On the other hand, the reason why, in spite of the lack of the resummation 
of quite important contributions, the SA turned out to be working well in the region B remained unclear until in 
Refs. 0, EE] we proved that the factors x a in the DGLAP-fits for the initial parton densities mimic the total 
resummation of the leading logarithms of x. Besides, in Ref. [3] we suggested to combine the DGLAP- results for gi 



with the results of Ref. [ll| in order to obtain an unified description of g\ in the regions A and B, without singular 
initial parton densities. A prescription for extending g\ into regions C and D was given in Ref. 161. 

In the present paper we present a unified description of gi valid in all of the regions A-D . The paper is organized 
as follows: in Sect. II we briefly remind the DGLAP-description of g\. For the sake of simplicity we consider in 
more detail, throughout the paper, the non-singlet component of g\ , and summarize the singlet results only. As 
the expressions for g± involve convolutions, they look simpler when an integral transform has been applied. The 
conventionally used transform is the Mellin one. However, in the small- a;, region, it is more convenient the use of 
the Sommerfeld- Watson transform, whose asymptotics partly coincides with the Mellin transform. This formalism is 
the content of Sect. III. Before dealing explicitly with gi, we consider in Sect. IV the appropriate treatment of the 
QCD coupling and compare it with the DGLAP-parametrization. The total resummation of the leading logarithms 
of x is quite essential in the small- a: region B. We discuss it in Sect. V by composing and solving appropriate 
Infrared Evolution Equations (IREE). Such equations involve new anomalous dimensions and coefficient functions 
and contain the total resummation of the leading logarithms of x. The singlet and non-singlet anomalous dimensions 
are calculated in Sect. VI. The non-singlet coefficient function is obtained in Sect. VII and is used to write down the 
explicit expression for the non-singlet gi in the region B. The singlet gi in the region B is obtained in Sect. VIII. In 
Sect. IX the small- a; asymptotics of the non-singlet g\ is discussed, whereas the singlet asymptotics is considered in 
Sect. X. Both asymptotic results are of the Rcggc type and their intercepts are found not so small. This may lead 
to the wrong conclusion that the IREE method cannot be applied safely to g\. In order to make this point clear, we 
discuss the applicability of our method in Sect. XI. In Sect. XII we compare our results for g\ in the region B to the 



1 Sudakov variables were introduced in Ref. 0] 
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DGLAP expressions. We show that DGLAP works well in the region B only because of the singular factors present 
in the parameterizations for the initial parton densities. On the other hand, when such fits are used, g DGLAP a j gQ 
behaves asymptotically as a sum of Reggeon contributions. The Regge behaviour in the two approaches is discussed 
in Sect. XIII. Combining the total resummation of the logarithms with the DGLAP results, we give in Sec. XIV the 
interpolation expressions describing g\ in the unified region A®B. Furthermore we show in Sects. XV and XVI how 
it is possible to describe g\ in the small-Q 2 regions C and D and arrive thereby to the interpolation expressions for 
gi which can be used in the whole region A © B © C © D. In particular the small-Q 2 regions C and D are described 
by a shift of Q 2 . Such shift is a source of new power Q 2 - corrections and we discuss them in Sect. XVII. Due to the 
experimental investigation of the singlet g\ by the COMPASS collaboration, in Sect. XVIII we give an interpretation 
to the recent COMPASS data. Finally Sect. XIX contains our concluding remarks. 



II. DGLAP -EXPRESSIONS FOR gi 

The Standard Approach to <?i is based on the DGLAP evolution equations and also involves some standard param- 
eterizations for the initial parton densities Sq and Sg. As the notations for the anomalous dimensions, the coefficient 
functions and the fits for the parton densities vary widely in the literature, we explain below the notation we use 
trough the present paper. 

We will denote gf s DGLAP and gf DGLAP the non-singlet and singlet parts of g\ when the SA is invoked. As the 
expressions for g\ involve convolutions, it is convenient to write them down in the Mellin integral form. In particular, 
the non-singlet g\ is: 



91 
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" l " :LAP {x,Q 2 ) = (e 2 /2) I £l(l/xrC» 8 DGLAP (uo,a s (Q 2 ))h{u>) 
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where C NS dglap {uj, a s (Q 2 )) is the non-singlet coefficient function, r y NS dglap (oj, a s (Q 2 )) is the non-singlet anoma- 
lous dimension and Sq(uo) is the initial quark density in the Mellin (momentum) space. With the one-loop accuracy 
(NLO) (see e.g. Ref. [3]), the expression for C NS DGLAP is 



rtNS DGLAP _ riNS DGLAP , a s{Q 2 ) n NS DGLAP 

Z7T 



with 



r<NS DGLAP _ 1 s-<NS DGLAP _ n 
^LO " 1 ' '-'NLO " 



Similarly, with two-loop accuracy, 



r 1 1 



where 



„NS DGLAP 



l [0) (n) = C F 



n 2 2n 2n + l 2 \2 n(l + n) 
a s{Q 2 ) to), x , fa s (Q 2 ^ 2 



W (o)t \ , / a *w V a)f \ 
In V 2n I 



.n(l + n) 2 

We have used the standard notations Si, 2 in Eqs. (|19|21|) : 



(18) 



^S 1 (n) + S 2 (n)-S 2 (n)\. (19) 



(20) 



(21) 



j=n J=n 

A 1 J A 1 J 



3=1 " 3=1 

They are defined for integer n. Their generalization for arbitrary n is well-known: 



(22) 



Si(n) = C + yj(n - 1), S 2 (n - 1) = — + tf(n) , 



(23) 
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with C being the Euler constant. The standard fits for the initial parton densities include the normalization constants 
N q g , the power factors x~ a , with a > and more complicated structures; for example, 

Sq(x) = N q x~ a {l - xf{\ + 1X S ) = N q x~ a ip(x). (24) 

All parameters N q , a, (3, 7, 6 in Eq. ([24]) are fixed by fitting the experimental data at large x and Q 2 « 1 GeV 2 . 

The expressions for gf DGLAP are similar but more involved and we do not discuss them in detail in the present 
paper. For the sake of simplicity through the paper we use g\ non-singlet for illustration, when it is possible, but we 
present the final expressions for both the non-singlet and singlet explicitly. 



III. SOMMERFELD- WATSON TRANSFORM 



As it is well known, the DGLAP expressions for g\ involve convolutions and in our approach we use them too. 
The standard way is to use an appropriate integral transform. Traditionally, the SA uses the Mellin transform. We 
will proceed slight differently. Our goal is to obtain expressions for g\ at small x and will start by considering the 
spin-dependent forward Compton amplitude T M „ related to W^ m as follows: 

Wff n = ^92^ ( 25 ) 

where the symbol 3 means the discontinuity (imaginary part) of X^„ with respect to the invariant total energy 
s = (p + q) 2 ) of the Compton scattering. At large s, when hadron masses can be neglected, 

s w 2pq(l - x) = w(l - x) , (26) 

so s w w at small x. The amplitude can be parameterized similarly to Wf^ n : 

= iMne^Xp—lSpTxix^) + (s„ - Pp ^-)t 2 (x, Q 2 )] (27) 
pql \ q A I J 

so that 

gi = -UT! , 52 = -Ut 2 . (28) 

We call Ti y 2 the invariant amplitudes. Exploiting the factorization, can be represented as the convolution of the 
perturbative and non-perturbative contributions (cf. Eq. In particular, 

Ti=T q ®5q + T g ® 6g (29) 

where 5q and Sg are related to Sq and 5g through Eq. (|28|) . In the Born approximation (cf. Eq. Q), 



T Bom _ e 2 * T Born = q _ /gpx 

q w — Q z ,c 



From the mathematical point of view, Eq. ([2T)]) as well as the DGLAP equations and the expressions for gi 
are convolutions, so an appropriate integral transform can be used. On the other hand, the phcnomenological Regge 
theory (see e.g. [l7]]) states that in order to study accurately the scattering amplitudes at high energies, one should 
use the Sommerfeld- Watson (SW) transform 18] . The asymptotic form of the SW transform partly coincides with 
the Mellin transform and often this form is especially convenient to account for the logarithmic radiative corrections. 
The SW transform is actually related to the signature invariant amplitudes T^' defined, in the context of DIS, as 
follows: 

T (+) = i[T( S ,Q 2 )+T(- S ,Q 2 )] , T(-> = 1[T( S ,Q 2 )-T(- S ,Q 2 )] (31) 

so that 

T(s,Q 2 ) = T (+) +T ( - ) , T{-s,Q 2 ) = T {+) -T ( -\ (32) 



Let us demonstrate that the signature of the Compton invariant amplitude T\ in Eq. ([2T)l is negative. Using Eq. 
we can represent T± in Eq. ([2T)l as the sum of the signature amplitudes T{ ' . In order to satisfy the Bose statistics, 
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T^v should be invariant to the permutation of the incoming and outgoing photons in Eq. (|27|) . i.e. to the replacement 
combining /i ^ v and q ^ — q. On the other hand, in the limit of large s, where the SW transform makes sense, the 
proton spin remains unchanged under such replacement because S p « Pp/M whereas s ~ 2pq — > — s. It immediately 
allows one to conclude that the amplitude T± should not be present in Eq. (|27|) . Therefore, 



9l (x,Q 2 ) = -^T[-\ s ,Q 2 ) . (33) 

The Compton amplitudes with the positive signature contribute to the structure functions F\ and Fi describing 
the unpolarizcd DIS. The calculation of the non-singlet component of F± and the non-singlet g\ is quite similar, so in 
the present paper we consider gf s in detail and give the results for Ff s in Appendix A. In order to account for the 
logarithmic contributions it is convenient to use the asymptotic SW transform for amplitudes in the following 
form: 

T<±)= P ^(4)>( U )^)( W)S ) (34) 

where y = \n(Q 2 //j, 2 ) and £ are the signature factors: 

£ (±) = -{e- 1 ™ ± l]/2 w [1 ± 1 + incj}/2 . (35) 

The integration line in Eq. (|34p runs parallel to 3w and S should be larger than the rightmost singularity of F^(lj, y). 
Quite often in the literature 5 in Eq. (f3"4")> is dropped. In the phenomenological Regge theory, the mass scale /i in 
Eq. (jM)) should obey p 2 -C s, otherwise it is arbitrary. We are going to specify it later in the context of g\. The 
integration contour in Eq. (|34[) , which includes the line parallel to the imaginary w-axis, as stated above, must be 
closed up to the left. Then the contour includes all ^-singularities of F^ (u>, y). As Eq. (j34|) partly coincides with the 
standard Merlin transform, it is often addressed as the Mellin transform and we will do the same through the paper. 
Nevertheless, we will use the inverse transform to Eq. (I34[) in its proper form: 

( S /p 2 ,y) (36) 



2 f°° 

F ( - ± \uj > y) = — / dpe^" 9fT (±) 

JO 



where we have denoted p = ln(s//i 2 ) . Eqs. ([36| and (|C2|I are supposed to be used at large s (s ^> p 2 ) where the bulk 
of the integrals comes from the region of small u) (co -C 1). Obviously, Eq. l|3"6"l) does not coincide with the standard 
Mellin transform. Finally Eqs. (|25I27I31|) lead to 



1 r°° du 



=- / — (-1) wF^fav) ■ (37) 



IV. TREATMENT OF a 3 AT LARGE AND SMALL x 



The rigorous knowledge on a s is provided by the renormalization group equation (RGE). According to it, the total 
resummation of the leading radiative corrections to the Born value of a s leads to the well-known expression 

1 (38) 



Mn(-s/A 2 ) 



where A = Aqcd and b = (UN — 2n/)/127r, with N = 3 and n/ being the number of involved flavors. Eq. (|38|) 
is the asymptotic expression valid at |s| ^> A 2 . The value of s in Eq. I|38p is negative. Eq. is often addressed 
as the leading order expression for a s and is obtained with the total resummation of the leading, single-logarithmic 
contributions. Corrections to Eq. (|3"8]) are also available in the literature 2 but we will not use them in the present 
paper because in practice the accuracy of the total resummations of the other radiative corrections, usually accounted 
for with various evolution equations, never exceeds the single-logarithmic accuracy. The minus sign at s in Eq. (|38|) is 
related to the analyticity: a s (s) should be real at negative s, but when s is positive, a s (s) acquires an imaginary 



2 For recent progress in RGE see e.g. the review 19]. 
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part. Conventionally, a s (s) at positive s is understood as the value of a s on the upper side of the s -cut. Therefore, 
— s = sexp(— vk) and 

W-)-- 1 1 _l/ ln(.s/A 2 ) + z^ x 

as[S) ~ 6[ln( S /A 2 ) - ztt] 6Vl n 2 (s/A 2) +7r 2j- ldyj 

Expressions (|38l39p are perturbative and asymptotic. In order to be consistent with the applicability of the pertur- 
bative QCD, \x defined in Eqs. (|13I14[) should be large enough: 

H > A. (40) 

An alternative way is to modify Eq. (f3"5|) in order to be able to investigate a s at s < A 2 . For example, there is the so 
called Analytic Perturbation Theory (APT) suggested in Ref. [2(j. It is based on subtracting from Eq. ([38]) its pole 
contribution at s = —A 2 . In the vicinity of the pole 



a s (s) 



1 r A 2 1 



6 In ((A 2 + \s\ - A 2 )/A 2 ) 



b 



A 2 2 



+ 0(M-A 2 ). (41) 



The result of the subtraction is called the effective coupling and is used instead of a s . Such a coupling can be used at 
any value of s. The recent results in this approach can be found in Refs. 2l|. However, APT does not allow one to 
get rid of the cut-off fi when the Sudakov contributions of the higher-loop Feynman graphs are involved. So, in the 
present paper we do not follow this approach. 

Now let us discuss how a s is incorporated into the expressions for the amplitude A of the forward annihilation 
of the quark-antiquark pair into another pair. The generalization to the scattering of gluons can be obtained easily. 
We assume that the external quarks are almost on-shell, with virtualities ~ /i 2 , keeping /i 2 <C s. In the Born 
approximation, the amplitude Asom is given by the following expression (see Fig. [I]): 



A B = -Ana (7(coO M (~P 2 )^ M (P 1 ) M '(P 1 )^ H, (~P 2 ) = U (-P^HMPl) U '(Plh^ U '(--P2) M Born^ 
s S + It s ' 

where s = (pi +P2) 2 ■ In Eq. (|42f and through the paper we use the Feynman gauge for intermediate gluons. According 



-Pi 



-P2 



Pi 



k=P\+P2 



t 



Pi 



FIG. 1: The Born amplitude Asorn- 



to Appendix A, the quark color factor C (co ' } = C F = (N 2 —l)/2N for the t -channel color singlet and C (co/) = -1/2N 
for the vector (octet) representation. Through the paper we mostly discuss the color singlet amplitude. We address 
M Born (s) as the invariant amplitude for this process in the Born approximation: 



M 



»(*) = -47ra s C (coi) 



s + le 



(43) 



By definition, a s in the Born approximation is a constant. The radiative correction to M? orn can be divided into 
two groups: 

(i) The corrections contributing to a s . 

(ii) The other corrections. 

Leaving the corrections to (ii) for the next Sects, we consider now the effect of (i). They transform the fixed a s 
in Eqs. (|42l43p into the running coupling. It is possible to fix the argument of a s , using the arguments of Ref. [22| . 
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Incorporating the radiative corrections from (i) to Asom leads, in particular, to insert the quark bubbles into the 
(horizontal) propagator of the intermediate gluon. In the logarithmic approximation, each quark bubble brings the 
contribution ~ nf In s. The gluon logarithmic contributions, each ~ N, come from more involved graphs but eventually 
all contributions lead to the factor b — (UN — 2ri/)/(127r) which multiplies the overall logarithm. Obviously, the 
argument of this logarithm coincides with the argument of the logarithm from the fermion bubble contribution and 
it is s. The total resummation of the leading radiative corrections from group (i) converts the fixed a s of Eq. 
into the well-known expression of Eq. (|38p and therefore converts the Born invariant amplitude M Born of Eq. 
intoM(°): 



M {0 \s) = -4ira s (s)C {col) ■ 



(44) 



In order to apply the Mellin transform to M™, we allow for the shift s — > s — fi 2 , with s ^> /i 2 , in Eq. ([44]) . Then we 
can write 



M™(i 



~ £(4Y>> M , 

-too 2tti \fi z / 



with 



UJ 

where A(uS) corresponds to a s (s) in the u> -space : 



A(co) 



dpe 



(45) 



(46) 



(47) 



In Eq. (|4"T|) we have denoted r\ = ln(^ 2 /Ag C£l ). The first term in Eq. (|4T|) corresponds to the cut of the bare gluon 
propagator while the second term comes from the cut of a s (s). They have opposite signs because of the famous 



anti-screening in QCD, which is the basis of the asymptotic freedom for a s . In the literature m! ^ and Ff V> are 



,(0) 



/ 



often called Born amplitudes (and we also follow this tradition ) in spite of the fact that they include the leading 
radiative correction from group (i). The radiative corrections to Mi , i.e. the corrections from the group (ii), are 

often included by using evolution equations. Such equations involve convolutions of , so they look simpler in 
the uj -space. We discuss this in detail in the next Section and focus now on the parametrization of a s in the parton 
ladders. As the treatment of a s for the color singlet Mq and octet My amplitudes is the same, we will not specify 
the channel below. 

First we remind that the well-known result a s = a 8 (Q 2 ) in the DGLAP equations follows from the parametrization 



a s (kl) 



(48) 



in every rung of the ladder Feynman graphs, where the ladder (vertical) partons can be either quarks or gluons. 
The notation k± in Eq. (148)) stands for the transverse components of momenta k of the vertical partons (quarks and 
gluons). The theoretical grounds for this parametrization can be found in refs. [H, [24], 25 1. The analysis of the 
parametrization of a s directly for the DGLAP equations was discussed in details in Ref. 26 1. In Ref. [12| we had 
shown that the arguments of Ref. [26[ in favor of using the parametrization l|48p in DGLAP can be used at large x 
only. Later, in Ref. (27j we made a more detailed investigation on this issue and showed that the parametrization (|48|) 
is always an approximation regardless of value of x. As the matter of fact, a s (k 2 L ) should be replaced by the effective 
coupling a e /f given by the following expression: 



irb 

= a s {n 2 ) + \ 
irb 



arctan 



[)" arctan (w^7A^) 



\n(k 2 JpA 2 )J Vln^/A 2 
arctan (irba s (k\/ I3)**j — arctan (irba s (fi 2 ) 



(49) 



where the longitudinal Sudakov variable (3 is defined in Eq. (|15p . However when the starting point /i 2 of the Q 2 
-evolution obeys the strong inequality 



p 2 > A 2 e n w 23A 2 , 



(50) 
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oigff can be approximated by the much simpler expression: 

at" » a s (kl/P). (51) 

If additionally x is large, ay? ~ a s (kj_). For practical use, the inequality in Eq. (|50[) can be expressed in terms of 
the discrepancy R{p) defined as 

|(lM)arctan(Vln(/z 2 /A 2 )) - a s (^ 2 )| 
JZ(/*) = ^ ■ (52) 

A simple calculation shows that R(fJ>) rapidly grows when /i decreases, ranging, for example, from i?(/i) = 5% at 
I? = 2800A 2 to R(/jt) = 10% at /i 2 = 250A 2 and i?(/i) = 50% at p 2 = 8.7 A 2 . As the DGLAP starting point of the Q 2 
-evolution is typically chosen close to 1 GeV 2 , the latter example shows that the DGLAP parametrization Eq. I[48p 
has an error of 50 % at such low scale . This statement is true for all DIS structure functions. In order to derive 
Eq. (|49[) we consider now the parametrization of a s in the integral expressions for the DIS structure functions. Here 
we partly follow the approach of Ref. [2^|. To this aim we consider the forward Compton amplitude T(x, Q 2 ) related 
to the structure functions by Eq. (p?5[) . Obviously, this equation is true for all DIS structure functions, so we drop 
here the signature superscript in T as unessential. One can show (and in this paper we will do it in the context of 
DGLAP and our Infrared Evolution Equations) that T obeys the following Bethe-Salpeter equation: 

r<*. ^ . r— + ,JJ0 ?J ^ S M ({9 + *)».*..«.) 4.^^. (53) 

The integral term in Eq. ([55]) is shown in Fig. O The notation M((q + k) 2 , k 2 , Q 2 ) in Eq. ([55]) corresponds to the 




FIG. 2: The integral contribution in Eq. (|53[) . The w -cut is implied, though is not shown explicitly. 

blob in Fig. [21 Besides the amplitude T, it can also include a kernel (splitting functions). The inhomogeneous term 
rpBom j g rp j n ^_ ne g orn approximation. To be specific, we consider the case when the horizontal parton in Fig.[5]is the 
virtual gluon with momentum p — k whereas the vertical partons, with momentum k, can be either quarks or gluons. 
When they are quarks, k 2 + it should be replaced by k 2 — m 2 + le, with m q being the quark mass, but this shift does 
not play any role for our consideration below. The factor 2wk 2 _ in Eq. (|53|) appears as a result of the simplification of 
the spin structure of the ladder Feynman graph in Fig. [2] We use the Sudakov parametrization (fT5|) for momentum 
k of the vertical partons. In terms of the Sudakov variables a, /?, k±, 

k 2 = —wa/3 — k 2 ^, 2qk = wxa + w(3, 2pk = —wa. (54) 

where w = 2pq. It is convenient to introduce a new variable m 2 = (p — k) 2 instead of a. Therefore, 

Using Eqs. (|54I55[) . we can rewrite Eq. ([53]) in a simpler way: 

T(x , . + J- 2 /J ikl f £ ^p^'/lb + *)», Vnf + «i, Q>) - (56) 
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The integration over ft and k\ in Eq. (|56p runs over the region 

fto < ft < 1, m 2 < k\ < w. (57) 
The value of fto follows from the requirement of positivity of the invariant energy (q + k) 2 of the blob in Fig. [5] 

fto~x+ ^ 2 - (58) 

Let us notice that the m 2 -dependence in Eq. ([55)1 can be neglected to the leading logarithmic accuracy that we keep 
through this paper. 

A. Integration in Eq. ()56[) at fixed a s 



Let us consider first the calculation of Eq. (|56|) under the approximation of fixed a s . From the analysis of the ladder 
Feynman graphs (see e.g. the review |9|) one can see that DL contribution comes from the region of large k\ where 

- k 2 w k\ > k\ = (3m 2 . (59) 

It allows to neglect the dependence of M on ftm 2 in Eq. ([56|) . It is convenient to integrate Eq. ([56| over to 2 , using 
the Cauchy theorem. The singularities of the integrand are: the double pole from the vertical propagators 

ftm 2 + k 2 ± -ie = (60) 

and the simple pole from the horizontal gluon propagator 

m 2 + is = 0. (61) 

The integration contour can be equally closed up or down. Traditionally (see e.g. Ref. Q) the integration contour is 
closed down which involves taking the residue at the simple pole (|61j) . so we arrive at the following result: 



T(x,Q 2 )=T B °™ + ^± f *jl f 1 dft(l-ft)M(ft,Q 2 ,k 2 ± ). 



(62) 



Obviously, when x is large enough, one can change the upper limit of the integration over k\ for Q 2 . Similarly, 
fto « x. Extracting the Born factor 1/ ft from M we can write 

M(ft,kl) = (l/ft)PT (63) 

where P is a kernel. To specify it, let us provide T with the quark and gluon subscripts through the replacement T 
by T r (with r = q,g). It leads to specifying PT in Eq. (|6"3"|) : PT = P rr /T r /. At last, assuming that we have used 
the planar gauge allows to identify (1 — ft)P rr > with the standard LO DGLAP splitting functions. Differentiation of 
Eq. (f6"2"| with respect to the upper limit of the A:^) -integration (which is Q 2 at large x) leads to the standard integro- 
diffcrcntial DGLAP equations for the Compton amplitudes T q ,T g . 

B. Integration in Eq. (156[) with running a s 

When a s is running, it is also convenient to use the Cauchy theorem for integrating Eq. (|56p over m . The 
integration contour can again be closed down. However, the spectrum of singularities in the lower semi-plane now 
includes the pole (|61[) and the cut of a s running along the real axis: 

fi 2 < m 2 - it < +oo. (64) 

Therefore, instead of Eq. (|62p we arrive at the more complicated expression: 



1 



2tt 



T(x,Q 2 )=T Bo ™ + —I dk{ I dft(l-ft) -f^M(ft,Q 2 ,ki) + (65) 



) 



A* 2 J Pa 

™^M(ft,Q 2 , ftm 2 + k 2 ± )Za s (m 2 ) ^\ 
m A {ftm 1 + k^Y 



i'- 
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where we have not used the assumption of Eq. (|59|) . The second term in the rhs of Eq. (j65|) is the result of taking 
the residue in the pole (|61[) and the third terms corresponds to accounting for the cut Eq. ([55]) demonstrates 

explicitly that it is impossible to factorize a a , i.e. to integrate a s (m 2 ) over to 2 , without making the approximation 
of Eq. (|59p . When this approximation has been made, we immediately obtain 

1 f w dk 2 r 1 

T(x,Q 2 )=T B ° rn + — -± / d(3(l-(3)M(P,Q 2 ,k 2 ± )a e Jf, (66) 



27r jy? fc i j po 

with a e Jf given by Eq. (j49|) . Indeed, in this case the integral over m 2 in Eq. (|66|) is 
1 f k ±! dm 2 1 1 



J = 



arctan 



6J M 2 to 2 ln 2 (m 2 ) + 7r 2 tt& L Vln^/A^ Vln(fc 2 //3) 



(w^a 2 !") ~ &vct&n iuhs))] ■ (67) 



Obviously, the term it 2 in the integrand in Eq. (|67|l can be neglected when /x obeys Eq. (|50|) . It leads to the 
approximative expression of Eq. (f5Tjl for cfjt . The approximation /3m 2 <§; A: 2 ^ was also made in Ref. [26| for the 
integration Eq. (|56[) over to 2 in the case of running a s . However, the integration contour in that paper was closed up 
in order to take the residue of the double pole at k 2 — (3m 2 + k\ — 0, which contradicts the assumption of Eq. ([59]) 
made in Ref. [26|. Taking this residue automatically led to the wrong conclusion that a e Jf = a s (— k\j (3) regardless 
of the value of [i. This error was found and corrected in Ref. 



V. DESCRIPTION OF g x IN THE REGION B: TOTAL RESUMMATION OF THE LEADING 

LOGARITHMS 

The region B is defined in Eq. (| 10[) . As it includes small x and large Q 2 , both logs of 1/x and Q 2 are equally 
important in this region and should be summed up. The most important logarithmic contributions to g\ in region B 
are the double-logarithmic (DL) ones, i.e. the terms 

~ a™ m 2n - fe (l/x) ln^Q 2 //! 2 ), (68) 

with k = 0, 1, ..,11, so they should be accounted for in the first place. Then the sub-leading, single-logarithmic (SL), 
contributions can also be taken into account etc. Therefore, an appropriate evolution equation for 171 in region B 
should account for the evolution both with respect to x and Q 2 while DGLAP controls the Q 2 -evolution only and 
cannot sum up the logarithms of x. In addition, the running coupling effects in Eq. (|68j) should be taken into account. 
To this end, a special attention should be given to the parameterizations of a s in region B and we will use here 
the results of Sect. IV. In order to resum DL contributions we use the alternative method of the Infrared Evolution 
Equations (IREE) , first suggested by L.N. Lipatov (see Ref. [28j|). Then it was applied to the elastic scattering of 
quarks in Ref. [1^ and in Ref. [3(3], with the generalization to inelastic processes (radiative e + e~ annihilation). Since 
then the IREE method has been applied to various problems and a brief review of the applications can be found in 
Ref. [3lJ. It is convenient to compose IREE for the Compton amplitudes T^ - ' related to gi by Eq. (|33|) . 



A. The essence of the method 



As we have mentioned above, the DGLAP -ordering of Eq. (fT3|) makes impossible to collect all DL contributions, 
regardless of their arguments, to all powers in a s . In order to account for them, the ordering of Eq. (I13p should be 
changed as in Eq. (fT4l) . This leads to the infrared (IR) singularities emerging from the graphs with soft gluons. In 
order to regulate them an IR cut-off /i should be introduced and therefore the result of such calculation becomes 
fi -dependent. The fermion (quark) ladders contributing e.g. to the non-singlet components of the DIS structure 
functions do not need an IR cut-off as long as the quark masses are accounted for. But in order to treat them 
similarly to the graphs with soft gluons, one can choose \i ^> masses of involved quarks. After that the quark masses 
can be dropped and the only remaining mass scale is fi. Generally, the value of \x is arbitrary, with one important 
exception: in order to use the perturbative QCD, fi should obey Eq. PD|) . 

On one hand, such flexibility can be used to resum DL contributions through the use of evolution equations with 
respect to fi, which is the basis of our approach. 

On the other hand, after such a resummation has been done, we arrive to a result which depends on this indefinite 
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parameter 3 . Of course, the problem of fixing the IR cut-off is not a new one. It has been known since long time ago, 
appearing first in QED, where fi is replaced in the final expressions by a suitable mass or energy scales. However, in 
the context of QCD this problem becomes more involved. Indeed, besides regulating the IR divergencies, \i acts also 
as a border line between Perturbative and Non-Perturbative QCD. Of course such a border is totally artificial from 
the point of view of the physics of hadrons . 

Below we will discuss how we fix the value of \x for the non-singlet (where \i = 1 GeV approximately) and singlet ( 
fi = 5.5 GeV) components of the structure function g±. 

The last point deserving a discussion concerns the possible dependence of our results on the way of introducing the 
IR cut-off: basically, different ways can lead to different results as it was shown explicitly in Ref. [Hj]. However, such 
a discrepancy appears far beyond the leading logarithmic approximation (LLA), that we keep through this paper. We 
will discuss now the technical details of our approach. 

B. The IREE for T 1 ^ in DLA 

We will consider, from now on, the invariant amplitude t[ \ defined in Eqs. (|2"7|l and (151)) and related to the 
structure function 171 by Eq. (|33[) . To simplify our notations, we drop both the subscript and superscript at t[ ' 
and denote T = t[~\ When the cut-off fj, is used for calculating Feynman graph contributions to T, this amplitude 
acquires the additional dependence: 

T = T(w,Q 2 ,^). (69) 
This amplitude is in the left-hand side of the equation in Fig. [3) Beyond the Born approximation, T depends on its 




FIG. 3: The IREE for the Compton amplitude T (_) . 



arguments through their logarithms, so we can parameterize it in terms of logarithms: 

T = T(ln( W / M 2 ),ln(Q 2 /M 2 )) =T{p,y). (70) 

Therefore, 

^ dT _dT + dT 
dp? dp dy 

Eqs. ([70)) and (|7T|) are the left-hand sides of the IREE for T (see Fig. [3]) in the integral and differential form respectively. 

The right-hand side of the IREE includes, in the first place, the Born amplitude T^ rn given by Eq. ([50)) . It corresponds 
to the second term in Fig. [3) In general, the other terms of IREE are obtained by factorizing the DL contributions 



3 We remind that DGLAP is free of this problem due to the different ordering of Eq. 1131 . 
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of the softest partons. It is well-known that the DL contributions are technically obtained from the integration over 
both longitudinal and transverse momenta, each integration bringing a logarithm. Logarithmic contributions from 
the integration over each transverse momentum come from the kinematic regions where those momenta differ from 
each other: ki _l 3> kj ±. This implies that the set of virtual partons (quarks and gluons) always contains a parton 
with minimal transverse momentum. In other words, the transverse phase space can be represented as the sum of 
sub- regions Di, each of them contains the parton with a minimal transverse momentum. We address such a parton as 
the softest parton even if its energy is not small and denote k± its transverse momentum. The IR divergence arising 
by integrating dk±/k± can be regulated with the IR cut-off fi: 

k ± > (jl. (72) 

Of course, \x should obey Eq. (j4"0|) to guarantee applicability of the Pert. QCD. It was shown in Ref. 0] that in QED 
the DL contributions can be of two different kinds: 

(A) : the Sudakov DL contributions coming from soft non-ladder partons; 

(B) : non-Sudakov DL contributions calculated first in Ref. [||. They arise from ladder Feynman graphs. This 
classifications stands also for QCD. The DL contributions of the softest partons from groups (A) and (B) are 
factorized differently. 

Factorization of the softest gluon from group A: 

The DL contributions from the softest non-ladder gluons can be factorized by using the QCD -generalization [30l[33l .[33| 
of the Gribov factorization theorem (often called the Gribov bremsstrahlung theorem) of the soft photons obtained 
in Ref. [35| . According to it, the non-ladder gluon having the minimal k± and being polarized in the plane formed by 
the external momenta can be factorized, i.e. its propagator is attached to the external lines only. When the Feynman 
gauge is used, the softest gluon propagator connects all available pairs of the external lines. The integration over 
other transverse momenta have k± as the lowest integration limit. Such a factorization deals with k± only and does 
not involve longitudinal momenta. Obviously, there is no way to attach the softest gluon propagator to the external 
lines of amplitude T^~' with the DL accuracy 4 . 

Factorization of the softest partons from the group B: 

Both the DGLAP-ordering in Eq. (fTSJ) and the ordering in Eq. (fT4")) imply that one can always find a ladder (vertical) 
parton (quark or gluon) with minimal transverse momentum. However, there is a difference between the two cases: 
the softest parton in (fl3|) is always the lowest parton at the ladder whereas in (fT"4f the softest gluon can be anywhere 
in the ladder, from the bottom to the top. Therefore, it corresponds to the factorization of the lowest ladder 
rung in the DGLAP ordering (fTB")) . and the factorization of an arbitrary ladder rung under ([H)) . The latter option 
corresponds to the last term in Fig. [3l By definition, in both cases the integration over k± involves fj, as the lowest 
limit whereas integrations over other ki ± are [i -independent. 

Now we can compose the IREE for T in the integral form. The lhs is just T while the rhs consists of the Born 
contribution and the term obtained by using the factorization B, therefore we arrive at the following IREE 

T r (p,y)=T rBorn + i J ( 2*)4 (fc 2 W +} e) 2 T r'((<l + fc) 2 > Q 2 > k 2 )M r , r ((p - k)\ k 2 ) (73) 

where any of r,r' denotes q or g; the factor 2wk\ comes by simplifying the spin structure; The negative signature 
amplitudes M r / r of the 2^2 -forward scattering of partons correspond to the lower blobs in the rhs of Fig. [3] This 
will account for the total resummation of the leading logarithms as we'll see in detail in the next Sect. When the 
DGLAP-ordering lfl3|) is used instead of (fT"4|) . only the lowest ladder rung can be factorized and therefore M^J are 
in this case given by the DL part of the LO DGLAP splitting functions. In other words, we arrive in this case to 
Eq. ([53]). Applying the operator -\?djd\i 2 to Eq. {73]) it converts the lhs into Eq. ([TTJ) . On the other hand Eq. ([30]) 
shows that the Born contribution in the rhs of Eq. (|73|) vanishes under the differentiation. because it does not depend 
on [i. Using the SW transform as in Eq. (j34|) we rewrite Eq. (|73| in terms of the amplitudes F r , related to T r through 
Eq. {35): 



u>F r {uj,y) + 



dF r (uj,y) 
dy 



1 

8^ 



F T i(w,y)L T i r (w) 



(74) 



4 We will use this kind of factorization in the next Sect, for calculating the lower blob in Fig. [3] 
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where £ r y is related to M r i r by the transform Eq. ([34]) . The derivation of Eq. (|74|) from Eq. ([73]) is given in detail 
in Ref. [ll(. The general technique of simplifying the convolution in Eq. ([73")) is given in Appendix C. It is useful 
to rewrite Eq. (|74|) in terms of the flavor singlet, T s , and non-singlet, T NS , components of the Compton amplitude 



dF NS (u;,y) , 1 _ 



<9y v 8tt 2 

and 



w + Sj^HpT^y) + 3=2 (w,»)i OT (w), (76) 



<9y V 8tt 2 tov V « v 8tt 



It is also convenient to introduce the amplitudes related to as follows: 

Hik = Life- (77) 

F s , and F ws are related to the singlet and non-singlet components of 31 with Eq. (|3"7|) . Eqs. (|75l76p are written in 
the DGLAP-like form, with the derivative with respect to Q 2 , but actually they combine the evolution with respect 
to Q 2 and w. 

Let us consider how to incorporate the single-logarithmic corrections from group (ii) of Sect. IV into Eqs. 1|75|76|1 . 



C. Inclusion of single- logarithmic contributions into Eq. (|74|) 

Technically, the DL contributions appear from the integrals over the loop momenta fej of the following form: 

^f^kl), (78) 

where we have used notations p r for external momenta and presumed that j 7^ i. The function <p(p r , kj) is independent 
of ki, which follows from imposing the strong inequalities giving rise to the DL integration region: 

h x < k,j ; _l, ft <ft. (79) 

When, for example, linear terms in fcj are present in ip, one of the integrations in Eq. (|78[) does not give rise to a 
logarithm and as a consequence a single-logarithmic (SL) contributions appears. This takes place in the integration 
region where the strong inequalities 1(75]) do not apply. In particular, when the inequality fc, 1 < kj ± is not fulfilled 
there is not a single soft parton in this region and therefore the method we use cannot account for such contributions. 
On the other hand, replacing the DL inequality ft -C ft by the single-logarithmic one, ft < ft is not essential for the 
method and these SL contributions can be taken into account. This replacement converts Eqs. (|75|76[) into 



- F X S U'- </) + dFN l^ ,y) = 8^(1 A„„^.,„U')^ X'. ,.. ,si„ 



and 



cuF q s (cu, y) + «AJ = (i + A^wJM") if (u;, y) + (1 + A OT w)M<") ^/(w. f)' ! s:l 1 



5 F^' (lj v) 

uF g s {Lu,y) + = (1 + A^M") Fffay) + (1 + A OT w)/i„(w) 

with X rr i given by the following expressions: 

X qq = l/2, X qg = -l/2, \ gq = -2, X gg = -13/24 + n//(12iV), (82) 
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obtained with one-loop calculations in the planar gauge. The new amplitudes hn~ are defined similarly to Hik and 
account not only for the total resummation of the DL contributions but also contain the resummation of the SL 
contributions. Therefore Eqs. (|8QI8ip account for the resummation of the leading (DL) together with sub-leading (SL) 
contributions to the w -evolution and for the resummation of the leading (DL) contributions to the Q 2 -evolution. 
For this reason they can be considered as the generalization of the DGLAP equations to the region (B). In contrast 
to the well-developed technology of resummation of double logarithms, no regular methods for the total resummation 
of SL contributions are presently available in the literature. 



VI. IREE FOR THE AMPLITUDES H ik AND h ik 



The expressions for the parton amplitudes Hik and hik in Eqs. (|75H81[) can be found by explicitly solving the IREE 
for them. Such equations were obtained and discussed in detail first in Ref. 1 101 ] where a s was kept fixed, while 
the running coupling effects were implemented in Ref. pT| . The technique for solving the IREE for the anomalous 
dimensions is similar to the one for amplitudes T , so we give below just a short comment on it. As it was done 
before for the Compton amplitudes, we begin with the IREE for amplitudes Mik(p) related to the Mellin amplitudes 
Hik(u}) through Eq. (j34|) . The amplitudes Lik(p) do not depend on Q 2 , so the left-hand side of the IREE for Hik{u>) 
does not involve derivatives and is equal to ojHik{to). The invariant amplitudes M in the Born approximation were 
introduced in Eq. (|44j) . It was shown that in the u> -space all of them are ~ 1/w, so we can write 



TrBorn 



a*ik/u, 



(83) 



with 



A(uj)C t 

27T 



A'{oj)Cf 



'<l!l 



"91 



2tt 



l 99 



WA(u) 
2^ 



(84) 



where we have used the standard notations Cf = {N 2 — 1)/2N = 4/3, rif is the number of the quark flavors, A is 
defined in Eq. (J4TJ) and 



A'M = - 



dpe~ 



iO p 



(p + v)' 2 



(85) 



The reason for the replacement of A(u>) by A'(u>) is that the argument of a s in L qg and L gq is space-like, so the coupling 
does not lead to tt -terms. The next term in the rhs of the IREE are the convolutions of the anomalous dimensions 
appearing as the convolution in Eqs. (|75I76[) . At last, the rhs contains the contribution from the factorization of the 
soft (Sudakov) gluons. We account for this contribution approximately (see Ref. [TTj] ) for more details. Eventually 
we arrive at the system of IREE for the singlet anomalous dimensions Hik represented in Fig. [H In the case of the 
non-singlet anomalous dimension, all gluon contributions in Eq. (|8"6")) should be dropped. Rewriting Fig.0]in a detailed 
form, we arrive at the system of the algebraic non-linear equations for Hik ■ 



uH qq — b qq 

luH, 



91 



J gq 



HqqHqq 
HgqHqq 



HqgHgq^ 
Hgg Hgq , 



UlH qg — bqg 
LOH, 



99 



J gg 



HqqHqg + HqgH gg , 

HgqHqg + HggHgg. 



(86) 



where we have denoted 



(87) 



with dik given by Eq. 



Then 



7T Z 



m qq = 2fi> m 99= ) TO 9<Z = > m l9 = ~NC F , 



(89) 



and 



1 f°° 

D (") = W J dpe-^lndp + ri/rj) 



P + V 



1 



(p + r\) 2 + 7T 2 p + 7] 



(90) 
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FIG. 4: The TREE for singlet amplitudes Hit . The solid lines correspond to quarks and the wavy lines denote gluons. 



Let us add a comment on the terms Vik ■ They appear from those graphs in Fig. 2] where the softest virtual gluon 
is factorized, i.e. its propagator is attached to the external lines in all possible ways. The blob obtained after the 
factorization is not a color singlet but an octet in the t (vertical) channel because the factorized gluon bears the color 
and belongs to an octet representation of the color group SU(3). So, these new amplitudes M\ k should be calculated 
independently. It is not a big deal for the non-singlet gi which involves M\ q only, but becomes a serious technical 

problem when all are involved (see Ref. Fortunately, all M^f 1 rapidly decrease with energy, and then is 

possible to approximate them by their Born values with a few per cent accuracy as was suggested in Ref. [3r| . 

Eqs. (|86p . similarly to Eqs. (|75l76p . combine the total resummation of DL contributions and the running coupling 
effects but do not include other SL contributions. The part of SL contributions accounted through Eq. (|82p can be 
easily incorporated into Eqs. (|86|) . leading to the following equations: 



Ujflq 
Luh q 
Ujh n 



b qq + (1 + \qqljJ)hqqhqq + (1 + \ qg Ul)h qg h 
b qg + (1 + \q g L))hqqhqg + (1 + Xqg 
bgq + (1 + X qq Uj)hgqhqq + (1 + \ qq Uj)hggh 



(91) 



99^99' 



Uhgg = b gg + (1 + \qqUj)hgqhqg + (1 + XggUj)h gg h gg . 

Now we present the expressions for the non-singlet case. Dropping all gluon contributions in Eq. (|91|) , we immedi- 
ately arrive to the equation for the non-singlet amplitude h qq : 



wh a 



bqq + (1 + XqqUj)(hqq) 2 



(92) 



with the obvious solution 



where 



[u> - \/u) 2 - B NS ] 
2(1 + \ qq uj) 



(93) 



B N S = 4(1 + X qq Uj)bqq. 



(94) 



Unfortunately, the system of non-linear algebraic equations in Eq. (f9Tj) can be solved analytically only if all Xik are 
dropped. In this case the solution to the system is 



H„ 



1 


U) 


-z + 


bqq 




b 99 


2 




Z 




1 




-z- 


bqq 




b 99 


2 




Z 





H n 



"<I9 

z 

Z 



(95) 
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where 



Z=^= ] J(^- 2(b qq + b gg )) + ^{^- 2(b qq + b gg )Y - A{b qq - b gg y - Wb gq b qg . (96) 

The non-linear algebraic equations in Eq. ([92]) and Eq. ([86]) have more than one solution, however the solution chosen 
in Eqs. (|93I95|) obeys the matching condition 

h m -> a qq /u, H ik -» Hg° rn = a lk /uj (97) 

at w — > oo, with given by Eq. (|84p . In other words, this matching condition in Eq. (|97p implies that these 
amplitudes are represented at low energies by their Born values. Indeed, from Eq. (|34[) high energies correspond to 
small u and vice versa. 



VII. SOLUTION TO THE IREE (SO) FOR gi NON-SINGLET 

As soon as the expressions for h qq are obtained, one can easily find the general solution to the linear differential 
equation ([80]) for the Mcllin amplitude F NS . As the procedure between the non-singlet and singlet cases has a purely 
technical difference, we consider in detail the former case and proceed to the singlet case in a much shorter way. 



A. General solution to the non-singlet equation (|80J) 

Obviously, the general solution to Eq. (|50|) is 

pNS _ pNS^ e -uy+y(l+\u)h qg (gg\ 



and therefore 



T NS (x,Q 2 ) = I' —( w /A UJ F NS {uj)e- bj y + y^ +x ^ h - = f°° —{l/ x YF NS {u:)e^ 1+x ^^ (99) 

J_ 400 2l7T V J J_ loo 2l7T 

with F NS (co) being arbitrary. In order to specify it, we use the matching condition 

T NS (x, Q 2 ) = f NS (w/n 2 ) (100) 

when y = 0. The new amplitude T NS (w/ /j, 2 ) describes again the forward Compton scattering off the same quark, 
however the virtual photon has now a virtuality » /j 2 . The IREE for T N s (w / /j, 2 ) should be obtained independently. 



B. Composing the IREE for T NS (w/fi 2 ). 

The IREE for f NS (w/fi 2 ) is similar to the IREE for T NS (x, Q 2 ), Eq. ([73]). It has the same structure and involves 
the same amplitude M qq . Still, it differs from Eq. (|T3"|) because of two following points: 

(a) T NS (w/ p 2 ) does not depend on Q 2 , so the differential IREE for it does not involve d/dy; 

(b) the Born amplitude Tg^ rn (w / /j 2 ) can be obtained from T NS (x, Q 2 ), putting Q 2 = fj,, so its contribution does not 
vanish under differentiation with respect to fi. Then introducing the Mellin amplitude F ns (lu) related to T NS (w//i 2 ) 
through the transform ([34]), we arrive at the following IREE: 

u;F NS (u,) = (e 2 /2) + (1 + \ qq uo)h qq F N s (uo) . (101) 

and therefore 



(102) 
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C. Expression for g^ s in the Region B 



Combining Eqs. (|102[ [98)) and p7|) and convoluting the perturbative expression with the initial quark density 
immediately leads to 



9? S 



du> / w \ w 



U)5q(uj) 



2m\^J u _ (1 + X qq u)h qq (w) 
u>5q(u) 



zoo r~ — y±- T s\qq\v )iv q q\ 

2 rlo ° duo 



-uy+yh q g 



(103) 



2 •/-,«, 2« Vx/ w - (1 + \ m u)h qq {u)) 
where 6q(u>) is the initial quark density in u> -space. Confronting Eq. (|103p to Eq. (fTTf it is clear that 



tiNs(u) = (1 + X qq uj)h qq (uj) = (1/2) w - \/uj 2 - B NS (u>) 



(104) 



is the new non-singlet anomalous dimension. It contains the total resummation of DL contributions together with 
the running a s effects and a part of SL contributions as explained in the previous Sect. Similarly, 



C 



2w 



NS 



(1 + \ qq w)fl qq (w) U-h NS (uj) LO + y^LO 2 - B ns (lu) 



(105) 



is the new non-singlet coefficient function. It is expressed through the anomalous dimension and therefore incorporates 
the same kind of logarithmic contributions. Eventually we arrive at the final expression for the region B of 

large Q 2 and small x: 



3 2 />zoo 
,NS(„ /o2\ _ 1 



2m 



(106) 



VIII. SOLUTION TO THE IREE ([76]) FOR THE SINGLET g L 

Eq. ([76| for gf can be solved in a similar way as the IREE for g^ s . First a general solution should be obtained 
and then constrained with a boundary condition. The general solution is easy to obtain: 



F q = e 



-uy 



c {+) e n {+)y + cr(-) e fi(_)» 



(107) 



2H, 



2H, 



<w 



where 



and 



x — H gg ~ H qq 



R - {Hgg - H qq ) 2 + 4:H qg H g 



(108) 



(109) 



We remind that the anomalous dimensions Hik are found in Eq. (|93|) . The exponents fi(±) are also expressed in terms 
of H ik ■ 



H qq + Hgg i X / (Hgg - H an ) 2 + AH an H, 



Finally the quantities and ' have to be specified. 

In order to constrain Eq. (|107f> we use the matching condition at Q 2 = fi 2 : 

F q (uj, Q 2 = /i 2 ) = C<+> + = F q (uj) , 

Fg^Q 2 = (i 2 ) = + C(->^-^ = F g {u>) 



(110) 



(111) 



2H, 



'HI 



2H, 



•hi 
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The new Mellin amplitudes F q . g (ui) correspond to the forward Compton scattering, when the photon virtuality is /j, 2 . 
They should be found independently. We again proceed by using new IREE for them. They have a structure similar 
to Eq. (|76|) . The difference is that the new equations do not contain derivatives with respect to y and account for the 
Born contributions 

pBom = < e q > pBom = q (112) 

where < e 2 > is the standard notation for the averaged e 2 . So, the IREE for amplitudes F qg are 



coF q s =< e\ > h qq {u) F^lu) + h n (u) F g s (u,y) , (113) 
coF g s = h qg (u) F q s (u) + h gg (u) fjf( W ) 



and the solution to Eq. (|113[) is 



Ff = < e 2 > - ^ , (114) 

q UJ 2 — U)(H qq + H gg ) + {H qq H gg — H qg H gq ) 

F n s = <el> Hgq 



U 2 - U}(H qq + H gg ) + (H qq Hgg - H qg H gq ) 



'<?<? ' "ggJ ~ y+iqq^^gg ±± qg ±± gq 
Combining Eqs. (|1 14[) and (|111[) . we obtain the explicit expressions for 

C (+) = < e 2 > ZHggHgg - (X - VR)(m - H gg ) 

q 2Vi?k 2 - L0(Hqq + Hgg) + (H qq Hgg - HqgHgq)} ' 

C (-) = <e 2 > -2HggHqg + (X + VR)(U - Hgg) 

q 2Vflk 2 - U(H qq + Hgg) + (HqqHgg ~ HqgHgq)} " 

Introducing the initial quark and gluon densities Sq and Sg respectively and using Eq. (|37p . we finally arrive at the 
following expression for gf in region B: 

»ffeQ 2 ) = 5/'°° |^(j)"[(c<+> e «™» + C!-> e n <->») u «u>)+ (116) 

'clil^Hrf-.t^,,),,,,) 
ltl qg ZJrl qg 



where Sq(u>) and Sg(uj) are the initial quark and gluon densities in the u> -space. C^Xuj) and Cg ±] (uj) are the singlet 
coefficient functions calculated in LLA. The exponents ^(±) ar e expressed in terms of H^ in the same way as the 
DGLAP exponents are related to the DGLAP anomalous dimensions. So, we conclude that Hik are the anomalous 
dimensions for gf in LLA. The total resummation of DL contributions to the singlet g\ under the approximation of 
fixed a s was done in Ref. [l(|. This result was used in Ref. [37j where a s was running: a s — a s (Q 2 ) according to 
the renorm group concept. However, in Sect. IV we showed that this parametrization does not stand at the small-x 
region and should be changed by the parameterizations of Eqs. (|47l85p . The application of the expressions for g\ in 
Eqs. (|106lll6p and the study of their impact on the Bjorken sum rule can be found in Ref. [38j . 



IX. ASYMPTOTICS OF THE NON-SINGLET gi IN THE REGION B 

The expressions Eqs. (|106|116|) represent g\ in the region B. Before discussing them in detail, we consider first 
their asymptotics at fixed Q 2 ^> and x — > 0. Strictly speaking, such asymptotics can be obtained by applying 
the saddle-point method. When the small-x behaviour is proved to be of the Regge type (power-like), one can use a 
short cut by finding the position of the leading (rightmost) singularity in the ui -plane. Of course, such singularities 
can be different for g^ s and gf and should be found independently. In the present Sect, we consider the small -x 
asymptotics of the non-singlet g\ . 
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A. Asymptotic scaling 



Let us assume that the initial quark density Sq in Eq. (I106[) is non-singular in x at x — > 0, so it does not contribute 
to the small- x asymptotics. Then by applying the saddle-point method to Eq. (| 1 06[) one deals with the non-singlet 
coefficient function and the anomalous dimension only. In this case the stationary point is (see Appendix E for details) 



UJ 



VB^s l + (l- K ) 2 (y/2 + l//B^) 2 /(21n 2 , 



with 



dVB/dw\ w=u 



(117) 



(118) 



and y = ln(Q 2 / M 2 ), £ = VQ7(*V)- 

It immediately leads to the Regge asymptotics for the non-singlets: 



with 



IEvs — 



2(1 - K)y/B^\ V2 (»/2 + 



7rV2ln 3 / 2 ^ 



(119) 



(120) 



When in Eq. (|11T[) y <C 2/^Bns{^) (let us notice in advance that at w = Ajvs it means that y <C 150 /i 2 ), the 
value of wq does not depend on y at all. Therefore Eq. (|1 19[1 can be rewritten as: 



with 



The factor cns is 



5 f S (x,Q 2 ) ~ {e 2 q /2)5q{LO Q )c NS T NS {0 



T JV s(0=r o/2 /ln 3/2 e 



2(1 - K)/(7TV^S) 



1/2 



(121) 
(122) 

(123) 



and does not depend on y. Eq. (|121[) predicts the scaling behavior for the non-singlet structure functions: in the 
region Q 2 <C 150 fi 2 , with depending on one argument £ instead of x and Q 2 independently. Therefore in this 
region 



g? s ~ gf s = n Ars (w )^(^)(Q 2 / a; 2 M 2 ) W0/2 
at x — > 0, with w being the largest root of Eq. (|125p : 



(124) 



w 2 -Biv S = 0. (125) 

We call the result of Eq. (|124[) asymptotic scaling: g^ 5 asymptotically depends on one variable Q 2 /x 2 only, instead 
of two variables x and Q 2 . The DGLAP prediction for the asymptotics of g? s in Eq. (|132[) is quite different. Below 
we compare these results in detail 

According to the results of Ref. [l4| (see also Appendix E), in the opposite case, when Q 2 > 150 fj, 2 , the rightmost and 
non- vanishing at w — > oo stationary point is again given by Eq. (| 1 25[> but the pre-exponential factor IIjvs essentially 
depends on y, so the asymptotic scaling in this region holds for g^ s jy. 

Let us notice that the sign of the non-singlet asymptotic behaviour is positive when 5q(u>o) is positive and coincides 
with the sign of g^ s in the Born approximation. In other words, both the x and Q 2 -evolutions do not affect the sign 
of gf s . Now we focus on solving Eq. (|125|) . 
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B. Estimate for the non-singlet intercept 



In a more detailed form, equation Eq. (|125[) is 

" 2 -< 1 + ^)[(^)( 

dpe- up \n((p + T])/ri) 



n 



(2C F y 
V nb ) 



1 



ANC F 



p + 7] 



dpe~ 



(126) 



(p + r,f 



= 



where we have used the notation r\ — ln(^t 2 /A 2 ). Let us remind that in the case of fixed a s Eq. (|126[) is much simpler 
and can be solved analytically: 



uj 2 - 



2a s C F (2a s C F \ 2 1 



^Ia s U F y 



LU 2 4:NC F 



0, 



with the obvious solution ujq L given by the following expression [l( 



,DL 



{2a s C F /^l\l- 



N 2 - 1 



1/2- 



(127) 



(128) 



In contrast, Eq. (|126[) cannot be solved analytically. Besides, there is a big qualitative difference between the cases of 
fixed and running a s . Although the IR cut-off p is used for regulating the IR divergencies in both cases, Eq. (|127l) is 
free of any p -dependence whereas Eq. (|126[) is obviously p -dependent and therefore the solution loq also depends on 
p: ujq = ojq(p). The value of p is restricted by Eq. (|40|) only. As a result, we arrive at the solution to Eq. (|126[) in the 
form of the curve plotted in Fig. [5] Eq. (|126p shows that ujq depends on p through r/, therefore uiq depends on the 




. 1 



10 100 
p/A 



FIG. 5: Dependence of the intercept luq on infrared cutoff /i : 1- for Fi ; 2- for <?i ; 3- and 4- for F\ and respectively 
without account of 7r 2 -terms. The structure function is discussed in the Appendix B. 

ratio p/A and on rif. Besides the ij -dependence, loq is not sensitive to the value of A. The plot in Fig. [5] shows that 
the curve coq — LUo(p) rapidly grows at p < A, however this region contradicts Eq. (|40[) . so the perturbative expression 
(13^)1 for a s cannot be used at so small p. Both Eq. (|126|) and the plot in Fig. [S] are consistent in the region (|4"0|) only 
and should not be considered out of this region. Eq. (|126p has one maximum in region (I40[) : 

Ans = max [w ] = ^o(pns) = 0.42 (129) 

at 



p = p NS = Ae 2 3 w 10A. 



(130) 
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For the sake of simplicity we chose in Ref. n/ = 3 and A = 0.1 GeV. It would have been more realistic to choose 
A = 0.5 GeV, and Eq. (|130[) shows that such a change of A leads to multiply by a factor of 5 the values of fi for the 
singlet and non-singlet obtained in Ref. [ll| . Comparison of the curves 1 and 2 in Fig. [5] to the curves 3 and 4 shows 
that the important role played by the ir 2 -terms in a s (i.e. respecting the analyticity) for producing a maximum in 
the curves and 2. Furthermore in the vicinity of this maximum, the power expansion 

w o(m) — ^oi^NS) H j- (M-Mjvs)+2 ^2 (V-^ns) +••• (131) 

does not contain the linear term, so cJo(a*jvs) = Ajvs is much less dependent on /i than all other points on the curves 
1 and 2. This remarkable feature allows us to identify Ajvs as the best candidate 5 for the perturbative estimate of 
the genuine intercept of the non-singlet g\ . According to the prediction of the Regge approach, the genuine intercept 
should be a constant, with no other dependence. However, Eq. (|126p and its solution (|129[) account for the leading 
logarithmic contributions only and leave aside sub-leading perturbative contributions and possible non-perturbative 
ones, so it is hardly possible to identify (|129|) with the genuine intercept. Nevertheless, it turned out that our estimate 
(1129p is in a good agreement with the results of Ref. obtained by fitting all available experimental data. This 
leads to a very interesting conclusion: by some unknown reason all sub-leading and non-perturbative contributions to 
the non-singlet intercept happen to be either small or irrelevant at /i = (J-ns, so that the LLA prediction (|129[) proves 
to be a good estimate for the non-singlet intercept. Motivated by this result, we call /ins the Optimal non-singlet 
mass scale. However, it is worth stressing that this scale is an artefact of our approach and should disappear when 
non-perturbative contributions (also dependent on the same scale) would be accounted for. To conclude, let us notice 
that the Regge form of the small- x asymptotics of g± is the direct consequence of the total resummation of logarithms 
of x and cannot appear at fixed orders in a s . This small- a; asymptotics depends on Q 2 only through the factor 
(Q 2 // i2 ) A ™ s/2 (seeEq. (TT24J)). In particular it means that the intercept Ans has no dependence on Q 2 . On the other 
hand, the well-known DGLAP small- x asymptotics 



sf ~^[^M 1Wk (»2)]. m 

can also be obtained (see Appendix F for detail) with the saddle-point method providing the initial parton densities 
are not singular at x — > (In Sect. XII we consider the alternative case presently used in the Standard Approach for 
the analysis of experimental data at small x). The DGLAP -asymptotics (|132|) clearly does not exhibit the Regge 
behavior. The same is true for the case where the anomalous dimensions and coefficient functions are calculated in 
high but fixed orders in a s which would correspond to the NN..NLO DGLAP accuracy (see Appendix G for detail). 
In principle, one might think that a generalization of Eq. (|132p could lead to the Regge asymptotics, however with the 
intercept depending on Q 2 . We show now that there are no theoretical grounds for such a scenario. Indeed, it follows 
from Eq. (|F7[) that the Q 2 -dependence in Eq. (|132[) is the consequence of the use of the DGLAP -parametrization 
a s = a s {k 2 _) and the DGLAP -ordering (flU)) . As explained in detail in Appendix F, at small x this ordering should be 
changed by the ordering (TTJJ . Then the upper limit Q 2 in Eq. (fT3"|) in the small- X region should be modified to w (see 
Ref. |27j for detail). After that the DGLAP asymptotics will not depend on Q 2 but at the same time will not have a 
Regge-type form. The Regge asymptotics is achieved by accounting for the resummation of the leading logarithms of 
x. It exhibits an asymptotic behavior much steeper than the DGLAP result (|132[) , not only with respect to x but also 
with respect to Q 2 . The comparison of Eq. (Tg^) to Eq. (fl2^|) shows that gf s DGLAPj g NS _> q when x Q . The 
question however arises: how small should x be in order to allow our asymptotic expression Eq. (|124|) to represent 
g± s reliably? We answer this question below. 



C. Applicability region of the small-a; asymptotics 



The asymptotic expression (|124[) for g^ s is obviously much simpler than the integral representation (|106[) and also 
much easier to work with. However, it is valid for very small x only. In order to determine when Eq. (|124|) reliably 
represents Eq. (|106[) . let us study numerically the ratio 



?NS{ ,o2N| _ ffj vs (a;, Q 2 ) 
9i S {x, Q 2 ) 



5 We are grateful to P. Castorina for this very useful observation. 
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at fixed Q 2 and different values of x. The result is plotted in Fig. [6] According to it, g^ s is reliably represented 
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FIG. 6: Rate of the approach to asymptotics for different Q 2 : solid curve 1 for Q = 10/x , solid curve 2 for Q — lOOfi , 
dashed curve for Q = fj,. 

by its asymptotic expression g^ s (x,Q 2 ) at x < 10 -6 only. So, strictly speaking, Eq. (|124[) should not be used at 
available values of x. However, in the literature one can find that Regge type (~ x~ a ) fits of the experimental 
data are used at much larger values of x, and such fits are reported to work well. We suggest a simple explanation 
to this: the phenomenological parameterizations including the Regge type fits have nothing in common with the 
expression in Eq. (|124[) obtained with the saddle-point method. In order to use such parameterizations at relatively 
large values of x (at x 3> 10 -6 ), one can choose the exponents a in the fits greater than the genuine intercepts. An 
analysis of such Regge parameterizations can be found in Ref. To conclude, we also notice also that sometimes 
Regge parameterizations are used with intercepts depending on Q 2 : they have no theoretical ground and contradict 

Eq. HMD. 



X. SMALL-i ASYMPTOTICS OF THE SINGLET gi IN THE REGION B 

The small- x asymptotics of gf in the region B can be obtained quite similarly to the non-singlet case, by applying 
the saddle-point method to Eq. (|1 16[) . The singlet asymptotics also exhibits the Regge behavior: 

gf (x, Q 2 ) ~ (llxy°{Q 2 /tfY°' 2 [A{uj Q )5q{uj ) + B(cj )S 9 (uj o )} (134) 

where ujq is the stationary point and A, B include the asymptotic form of the coefficient functions. The position of 
the leading singularity corresponds the largest root of the equation 

(u 2 - 2(b qq + b gg )f - A(b qq - b gg f - 16b gq b qg = 0. (135) 

Similarly to the non-singlet case, the position of the singlet leading singularity depends on /i, with one maximum 

ee A s = 0.86 (136) 

achieved at (we choose again rif = 3) 

l-i = Us ~ Ae 4 « 55A (137) 

which gives £ts/A w 55 GeV when A = 0.1 GeV. By repeating the arguments given also in the previous Sect., we 
call As the singlet intercept and call /xs the Optimal singlet mass scale. The Optimal singlet and non-singlet mass 
scales are quite different. Our perturbative estimate (|136[) is also in a very good agreement with the result obtained 
in Ref. [40( by fitting the experimental data. Eq. (|134ll shows that the asymptotic scaling is also valid for the singlet 
g\: asymptotically gf depends on one argument Q 2 /x 2 only. In contrast to the case of the non-singlet asymptotics 
(I124p . the interplay between 5q and 5g can affect the sign of gf. Indeed, in the Born approximation gf > but it 
can be negative (positive) asymptotically depending on the sign of A5q + B5g in Eq. (I134|) . 
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XI. APPLICABILITY REGION OF THE IREE METHOD 



In this Sect, we discuss the region of applicability of our approach and also answer a claim on a possible contradiction 
in our method that we have got in the past: in the IREE technology that we use to sum up the leading logarithms 
we work in the lu -space and systematically keep lu small. On the other hand when we calculate the non-singlet and 
especially the singlet intercepts, they are found not so small. Therefore we should guarantee the validity of our method 
not only at small but also at large u>. We start from the conventional analysis of the double-logarithmic QCD power 
series (|138[) . so we would like to stress at once that the use of (|138p for analysis of QCD processes and all estimates 
(e.g. the one in Eq. (|139[) ) based on it originate from the QED results (see e.g. Ref. [9|) where the running coupling 
effects can be neglected. They can become unreliable in QCD and therefore these conventional estimates should be 
replaced by more accurate estimates which we present in this Sect. 

Now let us remind the basic principle of the Leading Logarithmic Approximation, and DLA in particular. The 
straightforward calculation of the Feynman graphs contributing to a certain quantity (for instance, to the non- 
singlet coefficient) yields the double-logarithmic contributions. In the ro-th order of the perturbation theory they are 
~ a™ ln 2n (l/a:). For the sake of simplicity we keep here a s fixed and leave out other numerical parameters like 1/V, 
color factors, etc. Accounting for the running coupling effects and sub-leading logarithmic contributions does not 
change the essence of the problem. The series of such contributions, for example 

c\a s ln(l/x) + C2« s ln 3 (l/a;) + c 3 a^ln 5 (l/a;) + ... (138) 

converges when a s ln 2 (l/:r) < 1 only. It brings us to a rough estimate for the lowest limit for x: 

x > x min = exp[-l/^/oQ. (139) 

On the other hand, it is interesting to know the result of the total resummation of the DL terms at really small 
x < x m i n and even at x — > 0. The reason is that the DL terms are not so large compared to other contributions at 
x > x m in but dominate at small x. However, the series in Eq. (|138p diverges at x < x m in and cannot be summed up 
in this region. The solution to this problem is well-known: in the first place the series Eq. (|138[) should be summed 
up at x > x m in and then the result of the resummation can be analytically continued into the region x < x m in- Let 
us notice that the series Eq. (I138j) becomes divergent in region x < x m in and is called an asymptotic series. In the u) 
-space the series Eq. (|138l) , according to the relation 

°^e^) 1 1^(1/,), (140) 



is given by: 



This series converges when, roughly, 



2m uj 1+2n (2n) 



a _ a 2 „ a 3 



ci—+c 2 —+c 3 — + ... (141) 



U > Umin = \fa~s (142) 



whereas the DL terms becomes large at u> < u> m in ■ To be specific, let us notice that the expressions for the coefficient 
functions in Eqs. (|105I115|) represent the total sum of the DL contributions. Strictly speaking, the coefficient functions 
should first be calculated for large lu: lu > Lu m in (where the DL contributions are small) and then continued to the 
region of small lu. However, anticipating the analytical continuation into the small -lu region, quite often a short 
cut is taken and we follow this way: we treat lu as small since the beginning. It gives us the reason to neglect all 
non-logarithmic corrections regardless of their relatively large values in the region (|142p . After the total resummation 
of the DL terms has been done, the formulae obtained are insensitive to the value of lu and therefore can be used at 
any x. The resummed expressions (|105lll5p contain new singularities luq (branching points in our case but, generally, 
they can also be poles), which are absent in the series (|14ip . The rightmost singularities, i.e. the intercepts (Ajys and 
As in our case), determine the range of convergence of the series Eq. (|138p instead of ^Ja^: the series (|14ip converges 
only if 

lu > A, (143) 

with A being the intercept. However, after the total resummation in (|14ip has been performed, the result of the 
resummation can be used at arbitrary values of lu. To conclude, we note that all equations for the resummation of 
the leading logarithms, and in particular the IREE we have used, are not the equations for finding the intercepts. 
Indeed, the intercepts are the singularities and the values of lu in the IREE should be kept pretty far away from them 
by definition. The intercepts appear in the asymptotic expressions and therefore they should be found independently 
of the resummation methods, usually by applying the saddle-point method. 
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XII. COMPARISON OF gi TO pf IN THE REGION B 



In this Sect, we compare our results (|106I116[) to the DGLAP expressions for g\ at small x. We are not going to 
use here the asymptotic expressions in Eqs. (|124ll34p or (|132p , but we compare the two approaches at small but finite 
x. First we will compare the basic ingredients of the expressions for g\. the anomalous dimensions and the coefficient 
functions. Whenever it is possible, we will consider in detail, for the sake of simplicity, the non-singlet gi and more 
briefly generalize our results to the case of gf . 



A. Comparison of the coefficient functions and the anomalous dimensions 



Eqs. (|106p and (fT7|) have a similar structure: each integrand contains the initial parton density, the coefficient 
function and the exponent with the non-singlet anomalous dimension to govern the Q 2 -evolution. However, Cns 
and /ijvs in Eq. (|106p contain the total resummation of the leading logarithms of x whereas in Eq. (|17p the coefficient 
function and the anomalous dimension are considered to LO and NLO accuracy, namely they are given in Eqs. (|18l20p . 
Originally DGLAP was suggested for studying the region A of large x and large Q 2 . Due to the oscillating factor 
x~ u in the Mellin integrals, the main contribution to g\ in the region B comes from small u>. On the contrary, the 
main contribution in region A comes from large to. At large u>, the expressions for Cns and Iins in Eq. (|106|) can be 
expanded into a converging series in l/u>: 



Cns = 1 + 



A(w)C f 



2tt 



1 



1 

2uj 



(144) 



h-NS 



A(u)C F 
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(145) 



Obviously we observe a large discrepancy between Eqs. (|144ll45p and the LO DGLAP expressions in Eqs. (|18l20p . 
However, this discrepancy almost disappears when we come back to the region B where uj is small and therefore regular 
in Eqs. (|18l20p can be dropped. The remaining discrepancy is due to the different treatment of the QCD 



terms 



u> 



coupling. When the starting point of the Q 2 -evolution obeys Eq. (|40[) . then A(uj) with very good approximation 
can be replaced by a s (k 2 _/x), but definitely not by a s (k 2 _^). Taking into account more terms in the series and adding 
them to Eqs. (|144ll45p does not change the situation. So, we conclude that in region B the first and second terms of 
the 1/ui -expansion of Eqs. (|104ll05p reproduce the most important LO and NLO DGLAP results in the non-singlet 
anomalous dimension and coefficient function, with the exception of the different treatment of the QCD coupling. 
Expanding Eqs. (|95I115[) into a series in 1/ui and comparing the result to the singlet DGLAP anomalous dimensions 
and coefficient functions, we arrive at the same conclusion. 



B. Numerical comparison of the x -evolutions in Eqs. (|17p and (|106p 



The integrands in Eqs. (TT7|) and (| 1 06[) for the non-singlet g\ contains also a phenomenological ingredient: the initial 
quark densities Sq. Let us introduce the ratio 

9i S 

R NS = NS DGLAP ( 146 ) 

9i 

and study its x -dependence at fixed Q 2 , for example at Q 2 = 10 GeV 2 . Obviously, this cannot be done until 5q is 
fixed. The choice 

8q(x) = N q S{l - x) (147) 

corresponds to approximate the initial hadron by a quark and to neglect all influence of the hadron structure. Of 
course, such a choice cannot be used for phenomenological applications, but it makes possible to compare the x 
-evolutions in Eqs. (JT7J) and (|106|) . The substitution of the bare quark input into Eqs. (TT7]) and (| 1 06[) leads to the x 
-dependence of Rns plotted in Fig. 7. This shows that the impact of the leading logarithms becomes quite sizable at 
xq « 10~ 2 . So, we arrive to a sort of puzzle: 

According to the different behaviour of the x -evolution in Eqs. (fT?) and (|106p . the DGLAP -description of g^ s 
should have failed for x < 10~ 2 , but phenomenologically it is well-known that DGLAP works well at x < 10~ 2 . The 
solution to this puzzle is given below. 
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FIG. 7: Rise of Rns of Eq. (|146|l at small x in case of bare quark input and for Q — 10 GeV 



C. The role of the initial parton densities 

In order to clarify the problem, let us consider in more detail a standard fit to the initial quark density, as in 
Eq. IpOljl : 

Sq(x) = N q x- a {l - xf(l + jx s ) = N q x~ a ip(x). 

with all parameters N q , a, /3, 7, 6 being positive. As the fit is defined at certain fixed values of x = xq and Q 2 = Qq, 
all its parameters depend on xq, Qq. We define N q as the normalization. As the term x~ a — > 00 when x — > 0, we 
call it the singular term, although the fit is introduced at large x. We call <p the regular part of the fit because ip — ► 1 
when x — > 0. Once transformed into the to -space, the fit becomes a sum of pole contributions: 



Sq(uj) =N q (to -a) 1 + ^2 ( m k(u + k - a) + ~/(lu + k + 5 - a)' 



(148) 



k=l 



where = (i(J3 — l)...(/3 — k + l)/kl. The first pole in Eq. (|148[) corresponds to the singular term x~ a in Eq. 
We call it the leading pole. The other, non-leading poles in Eq. (| 148|) originate from the interference between ip(x) 
and x _a . Substituting Eq. (|148p into the DGLAP expression (fT7|) . we see that the contribution of the leading pole, 



~ns dglap tQ g NS dglap ig ( we drop the NLQ contribution here) 



9i 



NS DGLAP ( x ,Q*)= e JL Nq (\ 



C 



NS DGLAP 



/ ln(Q 2 /A 2 ) x7(Q) 
W Vln(^ 2 /A 2 W 



(149) 



Substituting the other terms of Eq. (1148|) into Eq. (fT7|) it leads to a contribution quite similar to that in Eq. (|149[) , 
however with a — > afe = a — fc and a — — 1. Obviously, a > afe. Therefore g^" 5 DGLAP in Eq. (|149|) is really the 
leading contribution to g^ s DGLAP at small x and actually it represents the small- x asymptotics of g^ s DGLAP , 
Confronting Eq. (|149[) to the very well-known expression (|132[) for the DGLAP asymptotics, we see that they are 
totally different. The singular terms are also included into the DGLAP parametrization of the singlet parton densities. 
It leads to the steep growth of gi at small x and provides the reason for the agrement between the DGLAP -description 
of the structure functions and the experimental data. Therefore the DGLAP success st small x is related to the use of 
singular fits for the initial parton densities. This is the solution to the puzzle. Now let us discuss the most important 
consequences of this result. 

First, let us confront the asymptotics of Eqs. (|149p and (|124[) . We see that the x -dependence in these expressions is 
identical: both formulae exhibit the Regge (power-like) behavior. It allows us to conclude that the singular term x~ a 
in the standard DGLAP fits mimics the total resummation of the leading logarithms of x. Therefore, the singular 
factors can be dropped when the total resummation of the leading logarithms of x is accounted for. In order to show it 
explicitly, let us study numerically Rns, using the standard DGLAP fit of Eq. (|24|) . The results are plotted in Fig. [8] 
We can observe that Rns is pretty close to unity only when the fit is used in the expression for g^ s DGLAP ^ 
whereas only the regular part, if, of the fit is used in the resummed expression (|106p . All other options drive Rns far 
away from unity at small x. So, the resummation of the leading logarithms leads to simplify the standard fits. Fig. [5] 
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FIG. 8: Examples of the small x -behaviour of Rns for different singular terms in the initial quark distribution fits of Eq. (I24|l : 
a — (curve 1), a = 0.576 (curve 2), a = 0.36 (curve 3). All curves correspond to Q 2 = 10 GeV 2 , and involve only a regular 
part, ip, in the parametrization of the initial quark density used for g\ of Eq. (|106p . 



explicitly demonstrates that the singular fit (|2~4")l and the total resummation of the logarithms lead to close values 
of <7i in region B. From a practical point of view, the use of the resummation is preferable because it allows one to 
construct new fits with a reduced number of parameters. From the theoretical point of view, the resummation is even 
more preferable. Indeed, the DGLAP intercept a in Eq. (|149p depends on the starting point Xq, Qq, where a is fixed, 
and such a dependence can hardly be deduced from theoretical considerations. On the contrary, the intercept Ajys 
in Eq. (|124[) is independent of the initial parton densities. 

Now we would like to comment on an apparent puzzle arising first from the results of Refs. [H, l4l| and then also in 
other subsequent publications (see e.g. the recent review [Hj]): 

On one hand, the direct comparison in Ref. [4l[ of the NLO DGLAP result, g^ s DGLAP to the expression for g^ s 
obtained in Ref. [l(| in the limit of fixed a s showed that 

g NS /g NS DGLAP <<; j (15Q) 

at the small values of x available in experiment and therefore the conclusion was made that the resummation of 
ln fe (l/x) can yield a small impact. On the other hand, it is clear that the small- x asymptotics (|124[) is much steeper 
than the well-known DGLAP asymptotics (|132[) . which proves that asymptotically 

g NS /g NS DGLAP ^ j (m) 

and therefore the total resummation of ln fc (l/a;) should be essential. So, Eqs. (|150[) and (|15ip obviously contradict 
each other, which is puzzling. 

Eq. (|150p was interpreted in the literature as follows: the resummation of leading logarithms at values of x available 
in experiment is much less important than the impact of the sub- leading (compared to the double-logarithmic 
contributions) terms in the DGLAP coefficient functions and anomalous dimensions. 

However, a close inspection of Eqs. (|24I148|) for the standard fit suggests another solution to this puzzle. Fig. 8 
clearly demonstrates that the main impact on the small- a; behavior of g^ s DGLAP comes not from the NLO DGLAP 
coefficient functions but from the singular term x~ a in Eq. (|24p . Indeed, when this factor is removed from Eq. (124p . 
we arrive at curve 1 despite the sub-leading contributions are accounted in g" 8 DGLAP . Therefore, their impact leads 
to Eq. (|15ip instead of Eq. (| 1 50[) . On the contrary, when the singular term is accounted for, we arrive at curve 2 
where Rns ~ 1- It proves that the conclusion of the extreme importance of the sub- leading contributions on the 
small- a; behavior of g\ advocated in Refs. [37],|4l[ is groundless. Now let us compare the small-x asymptotics of g± s 
and g^ s DGLAP . Parameters a in the DGLAP fits obey 



a > A, 



(152) 



28 



with A in Eq. (|152|) being either the non-singlet or singlet intercept, depending on the case. Eq. (|152[) naturally 
leads to Eq. (|150[) . When the singular factor x~ a in Eq. (fM| is dropped, the small- a; behavior of g P)GLAP is given 
by Eq. (|132p and the strong inequality sign in Eq. (|150j) should be reversed. The reason why the exponents a in 
the singular factors of the DGLAP fits should obey Eq. (|152p is clear: indeed, we have just shown above that the 
asymptotic regime is actually achieved at very small x, so in order to reproduce it at values of x accessible at present 
experiments, the parameter a, playing the role of the intercept, should be larg er than the intercepts A. In this 
connection we remind that our predictions agree very well with results of Refs |39tl4fH , whereas the phenomenological 
value of the intercept a in Eq. (|24p contradicts those results. It is clear that combining singular fits with the total 
resummation of logarithms also implies a double counting of the same logarithmic contributions: explicitly in the 
first case and implicitly in the latter, through the singular factors x~ a . Furthermore, Eqs. (|124[) and Q149p explicitly 
show that neither the DGLAP intercept a nor our intercept A^vs depend on Q 2 . Such a dependence, sometimes 
appearing in the literature as a possible generalization of Eq. (|132p is an ad hoc assumption and never appears as a 
result of QCD calculations. 

Finally, let is notice that it is commonly believed that the expression for the fit in Eq. (f24|) mimics the effect of the 
hadron structure, including basically unknown non-perturbative contributions. On the other hand, when the leading 
logarithms are accounted for and the initial parton densities are fitted at not too large x, the x -dependent terms in 
<p can be almost dropped, so the fit can be simplified down to N q . It means that the impact of the non-perturbative 
contributions is greater at large x whereas in the small- a; region it is reduced to a simple normalization. 



XIII. REGGEON STRUCTURE OF gi 

According to the Regge theory (see e.g. Ref. [l7:]), any forward scattering amplitude, including the invariant 
Compton amplitude T related to g\ through Eq. (j3"3"| . asymptotically exhibits the Regge (power-like) behavior and 
can be written as a sum of such power-like terms called Reggeons. The same should be true for g\. In this Sect, we 
show that both the standard approach (SA) and our description of g\ agree with such a representation. However, 
the Reggeons in these two approaches are different and the reasons for this Reggeon representation are also quite 
different. As usually we begin with considering in detail the non-singlet g±. 

A. Reggeon structure of gi in the SA description 
Eq. (|148p with the standard DGLAP fit in the ui -representation can be re-written as: 

oc oo — 

where r, r k , r k and j, j k , jk are expressed through the parameters of the fit as follows: 

j = a, jk = a - k, j k = a-k- 5, (154) 

r = N q , r fe = (l + 7)7V g /3(/3-l)..(/3-fe + l)/(fc!), r k = -yN 9 p(J3 - l)..(/3 - k + 1)/(W). (155) 

By inserting Eq. (|153[) into Eq. (fTT)) . integrating over u, and taking the residues of the poles of Eq. (|153[) . allows us 
to write g^ s DGLAP in the region B as the following series: 



NS DGLAP i 
HI 



where, to the LO accuracy, the terms S, Sf., St are 



2 00 

(a, Q 2 ) = f [S(x, Q 2 ) + ]T (S k (x, Q 2 ) + S k (x, Q 2 ))] (156) 



= (i) 3, c ™ 04 )^(^^)™°"' a ' ,/(2 " >) . 



29 



It is clear that the x -dependence of each of the terms S, Sk, Sk is Regge-like, so we call them the DGLAP Reggeons 
contributing to the non-singlet gi . The intercept j of the Reggeon S is the largest, and we call S the leading Reggeon 
and address Sk and Sk as the sub-leading Reggeons. Only the leading Reggeon has the positive intercept. All other 
intercepts are negative. We remind that all features of the DGLAP Reggeons are due to the assumed form of the 
initial quark density and are related to the phenomenological parameters of the fit (|24[) . Obviously, one can decompose 
the DGLAP expression for the singlet g\ quite similarly into a set of Reggeons. 



B. Reggeon structure of Eq. (|106|) 



Let us consider once more the limit of g\ at x — > 0. In Sects. X, XI we have shown that the use of the saddle-point 
method to Eqs. (|106I116[) led to the Regge asymptotics (|124ll34p . The intercepts Ams and Ag were determined in 
Eqs. (|129ll36p as the largest roots of Eqs. (]126[) and (|135[) respectively. They are not simple poles in the 10 -plane 
but the rightmost square-root branching points. They were found by solving numerically Eqs. (|126I135| . However, 
each of this equations can have more than one root. Applying the same argument we are able to find the additional 
non-singlet and singlet intercepts A^ s and Ag , with k = 1,2,... Accounting for these contributions allows us to 
represent the non-singlet and singlet g\ in a form similar to Eq. (| 1 56|) : 



9i {x,Q 2 ) 
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B{x,Q 2 )+Y,B k (x,Q 2 



(158) 



where the leading contribution B(x,Q 2 ) is given by Eq. (1124[) for g^ s and by Eq. (|134p for gf . The other Reggeons 
B k (x,Q 2 ) look quite similarly. Namely, they can be obtained from Eqs. (jl24[) and (| 13411 with the replacement 

Ajvs 



A$j and A s 



A§ respectively. In particular, for g^ b we have 



B k =U(AP s )Sq(A^ s )(Q 2 /x 2 ^) 



(159) 



and Reggeons for gf have the structure of Eq. (|134[) . We call B, B k QCD Reggeons because they are obtained from 
the total resummation of the leading logarithms in the QCD perturbation series. The Reggeon B has the maximal 
intercept compared to B k , so we call it the leading QCD Reggeon and Reggeons B k are the sub- leading (secondary) 
QCD Reggeons. It turns out that only the leading non-singlet Reggeon has the positive intercept (|129p whereas the 



next non-singlet intercept is A 
A^ = 0.55, A^ 2) = 0.35, a£° 



(i) 

NS 



0. On the contrary, there are three singlet Reggeons with positive intercepts: 



0.21. 



C. Comparison between the DGLAP and the QCD Reggeons 



The Regge theory, in the DIS context, states that the Regge (power- like) form of g\ should be achieved at x — > 
only, while g\ looks quite differently at large x. It perfectly agrees with the features of the QCD Reggeons B, Bk 
obtained with the saddle-point method from the expressions for g\ due to the QCD radiative corrections. They 
appear as a result of the total resummation of the QCD perturbation series and are never present to any fixed order 
of the perturbative expansions, including, of course, the Born term. Also they are not simple poles but square-root 
branching points, Their intercepts are found in terms of the basic QCD constants as the number of the colors N, the 
number of the flavors n / , and Aqc d ■ 

On the contrary, the SA Reggeons are produced by the poles present in any fixed order in a Sl including the Born 
approximation. They exist at any x, even at x ~ 1, because they are generated by the structure of the fit for Sq instead 
of QCD radiative corrections. The intercepts of the SA Reggeons are expressed in terms of the phenomenological 
parameters of the fit (f2~4"|) and have nothing to do with QCD calculations, so we call them input Reggeons in contrast 
to the QCD Reggeons B, B k . On one hand, the existence of such Reggeons contradicts the concepts of the Regge 
theory, On the other hand, we have shown that the phenomenological success of DGLAP at small x is due to the 
singular factors x~ a in the fits for the initial parton densities which mimic the total resummation of the QCD radiative 
corrections. Obviously, these parameters are chosen to match the experimental data. So one should not be surprised 
that a truncated set of input Reggeons S*, Sk, Sk could be close to the experiments with a good accuracy, the 
agreement being entirely due to the choice of their phenomenological parameters. So any theoretical interpretation 
of such Reggeons in the QCD context would be groundless. 
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XIV. DESCRIPTION OF gi IN THE UNIFIED REGION A0B 

In this Sect, we construct a description of g\ valid in both regions A and B. Again, we focus on the non-singlet 
gi in the first place. To begin with, let us remind that in region A , where x is large, the non-singlet gi is described 
by the DGLAP expression Eq. (JTTJ) where both the coefficient function and the anomalous dimension are known to 
the NLO (two-loop) accuracy. In order to describe the small- x region B, we took into account the leading 

logarithms of x and arrived at Eq. (|106[) . When Eq. (|106[) is considered in the region A , the LL contributions become 
small. On the other hand, non-logarithmic contributions accounted for in Eqs. (|18|20|) are quite important in this 
region. So, a possible option is to create an interpolation formula for g^ which would coincide with Eq. (|17[) and 
Eq. (|106|) in regions A and B respectively. To this aim, let us define new coefficient function C NS and anomalous 
dimension C NS by combining directly the DGLAP results of Eqs. (|18I20|) and the LL results of Eqs. (|105I104|) : 



C NS — C NS + C NS DGLAP 
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Because of the obvious double counting in Eq. (|160[) . let us make the necessary subtractions and define C^ ^ b and 
^comb' which we ca h the combined coefficient function and anomalous dimension' 
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where AC NS and Ah NS are the first- and second- loop terms of the expansion of C NS and h NS into the series (see 
Eqs. (114411451) ): 



AC NS = 1 
Ah 



A{uj)C f r 1 



2tt 

\ - _ A(uj)C f 



2n 



rl 

Ul 



1 

2^ 



(162) 



By inserting and in Eq. (|106[1 . we arrive at the final expression for g^ s valid in the region A0B 



(163) 



Quite similarly we obtain the combined coefficient functions C^^AuS) an( i anomalous dimensions h^ mb for the 
singlet gi : 



^comb ~ ° 
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ACW, h? k mb = h ik + hg GLAP - Ahi. 



(164) 



where C^q LAP correspond to the DGLAP coefficient functions with the replacement a s — ► A{lu) and hf k GLAP 
are the DGLAP anomalous dimensions with the same replacement. The subtraction terms in Eq. (I164p to the LO 
accuracy are: 
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Replacing and hik in Eq. (|1 by C$^L b and h^ mb , we finally obtain the expression for gf valid in both 

regions A and B: 
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where X and R are also expressed in terms of h^ mb . In the sub-region A of A0B the main contribution 

in Eqs. (|163I166[) comes from the DGLAP terms in the coefficient functions and anomalous dimensions while the 
logarithmic terms are small, so that Eqs. (|163ll66p almost coincide with the DGLAP expressions. On the contrary, 
in the sub-region B the main role is played by the LL terms and therefore Eqs. (|163I166[) are pretty close to the 
expressions of Eqs. (|106I116|) . Therefore Eqs. (| 1631 166|) really represent the interpolation expressions for g\ in region 
A0B. 

XV. DESCRIPTION OF g x IN THE REGION C 

The small Q 2 -region C is defined in Eq. [11] Contrary to the regions A and B, the SA cannot be used in region 
C at all. Indeed, the basic ingredient of SA, the DGLAP evolution equations, control the evolution with respect to 
ln(Q 2 /yU 2 ) in the regions A, B and do not apply at small Q 2 . In Ref.flfl] we have proposed a method to describe g\ 
at small Q 2 , which is a kinematic region studied experimentally. It turned out that our results for gi in region B can 
be generalized into the region C by introducing the shift 

Q 2 -> Q 2 = Q 2 + n 2 (167) 

where /j, is the infrared cut-off. Numerically, we have suggested to use the Optimal mass scales jj,ns = 1 GeV and 
(is = 5.5 GeV for the non-singlet and singlet case, respectively. The reasons for introducing those scales were given 
in Sect. X. Other shifts in Q 2 similar to Eq. (|167p were suggested in various papers, see e.g. Refs. [43|,|44j]. In the 
literature, such shifts were introduced from phenomenological considerations whereas we suggest it from the analysis 
of the Feynman graphs contributing to g\. Let us notice that introducing this shift we go beyond the logarithmic 
approximation we have kept so far, so in this sense we consider our description of gi in the region C model-dependent. 
To begin with, let us notice that both the singlet and non-singlet component of g\ obey the Bethe-Salpeter equation 
shown in Fig. [9l In the analytical form this equation is written 



q 




v 



FIG. 9: The Bethe-Salpeter equation for gi in the region C . 



where n stands for the numerical factors e 2 /2 and < e 2 /2 > for the non-singlet and the singlet respectively. We have 
skipped in Eq. (| 168[) the convolution with the initial parton densities in order to prove the shift in Eq. (|167[) . that 
the densities cannot affect. The S -function (together with the factor —2m) corresponds to the cut propagator of the 
upper quark with momentum k and mass m q coupled to the virtual photon lines and the factor 2k\ appears after 
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simplifying the spin structure of the equation. Similarly to Eq. (|33|) . E(2pk, fc 2 ) in Eq. (|168[) is related to the invariant 
amplitude M(2pk, fc 2 ): 

E(2pk, fc 2 ) = (l/27r)3?Af(2pfc, k 2 ). (169) 

The invariant amplitude M(2pk 1 k 2 ) describes the forward scattering of partons, with the upper partons being quarks. 
In other words, M(2pk, k 2 ) can be any of M NS (2pk, k 2 ) (for #f s ) or M qg (2pk, fc 2 ), M qq (2pk, k 2 ) (for the singlet gi). 
These amplitudes incorporate the total resummation of the leading logarithms. Obviously, E in Eq. (|168p does not 
depend on Q 2 . Contrary to the parton amplitudes entering in Eq. (|73[) . where k stands for the softest momenta, 
the amplitudes E in Eq. (|168p are essentially off-shell, and depend on two arguments: the invariant total energy 
(p + fc) 2 w 2pk and the virtuality fc 2 which is not small now. Therefore they should be calculated independently. The 
amplitudes M NS {2pk 1 k 2 ) 1 M qg (2pk 1 k 2 ), M qq (2pk,k 2 ) are considered in detail in Appendix C. 



A. Infrared regularization in Eq. (|168|) 

First, we remind that in order to account for the LL contributions in the small- x region C, one should use the 
ordering (fT4|) . This leads to the IR singularities of soft gluons and therefore an IR cut-off in the divergent propagators 
must be introduced. In order to treat the quarks and gluon ladders similarly, we drop the quark masses and introduce 
the IR cut-off in the quark and gluon ladders, providing both the ladder (vertical) partons and the soft non-ladder 
gluons with the fictitious mass /it, assuming that /i > m quar k- In the previous Sects, we introduced /i in somewhat 
different way: according to Eq. (|72[) . /i is the lowest limit in the integrations over k±. Although it was noticed in 
Ref. [32j that different ways of introducing the IR cut-off lead to different results, this goes well beyond the accuracy 
we keep. Let us also notice that there is no need to introduce /i into the horizontal propagators of the ladder because 
they are IR stable. Then Introducing /i Eq. (|168[) is modified into 

R „„ f d 4 k , ,,, N 2wk 2 E(2pk,k 2 + u 2 ) , s 

*=*? m +*J (2^ ( - 2 ^ ((g+fc)) (fc 2 -, 2 ) 2 2 Pfc ( 17 °) 



B. Solving the Bethe-Salpeter equation (|170[) 

It is convenient to write Eq. (|170|) in terms of the Sudakov variables defined in Eq. (fT5|) . Eq. (|D5|) shows that the 
2pk and fc 2 -dependence for any of E NS , E qq , E qg looks much alike, so below we consider the Bethe-Salpeter equation 
for g^ s only. Substituting E NS into Eq. (|170[) and changing the order of the integrations, we arrive at 
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In the region C x is small, so we drop the second term in the argument of the 8 -function, which is used for the 
integration over (3. We obtain 



gi = g?° rn + k 
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The region of integration in Eq. (|172p is shown in Fig. [TUl It is restricted by the following limits: 
(a): w > k\ + fi 2 > aQ 2 > 0; (b): wa > aQ 2 + k\ + \? . 

The result of the integration over this region depends on the relations between Q 2 and k\j 'a. The leading contribution 
comes from the sub-region D in Fig. 1101 Integrating over a in D leads to 



„ Born . , 

91 = fJi + « 



r°° duo uNS . . i r dk 2 ± i w \-/fc 2 +^\^ va 

; / — ojh NS (oj)- / , 2 ± „ 73 ? x , . (173) 

Replacing fc 2 + by t in Eq. (|173[) . we get 



= flf °™ + - / ^h Ns (c) / %wtnt/w HB . (174) 
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FIG. 10: The integration region in Eq. (|172j) . 



In the region Cbj>/i 2 ,so the upper limit of integration in Eq. (|1T4[) can be approximated by w. The lowest limit 
is definitely Q 2 + /i 2 which proves the validity of the shift we have suggested in Eq. (|167[) . Performing the integration 
over t in Eq. (|174p leads to the expression of Eq. (|106|) but with the shifted value of Q 2 . Indeed, the integration over 
t yields 

Born r°° du h NS r/«,NW/<32 N /.*s ,,„ x fc 



2m (tu - h NS ) 



The integration of the second term in the squared bracket yields zero, by closing the integration contour to the right 
of the singularity of h NS / (u> — h NS ). Using the identity 

(176) 



(u-h NS ) (uj-h NS ) 
and noticing that the first term in Eq. (|176[) cancels the term gf orn in Eqs. (|1681 IT75]) . we arrive at the following 



expression: 



2m (u - h NS ) \xJ V/i 2 . 

which coincides with <?f s from Eq. (|106[) with the replacement Q 2 — > Q 2 . We have used here the shifted variable 

x = Q 2 /w = x + fi 2 /w = x + z. (178) 

So, we have proved that our result for g^ s in region B can be extended to region C with the shift Q 2 — * Q 2 = Q 2 + . 

It is not difficult to repeat the above calculations for the singlet g\. Eventually we conclude that our expressions 
(|106I116|) for gi in region B (= g^\x, Q 2 )) can represent g\ in the region C with the shifts Q 2 — * Q 2 = Q 2 + fi 2 , x — > 
x, = x + z: 

g[ c \x,z,Q 2 ,f, 2 ) = g^ ) (x,Q 2 ). (179) 



Obviously, Eq. (|179|) is valid in the unified region B © C. 



XVI. DESCRIPTION OF gl IN THE FULL REGION A ffi B © C © D 

In this Sect we will show that combining the shift of Q 2 introduced in Eq. (|167|) and the interpolation expressions 
of Eqs. (|163II166|) for g\ in region A, it allows us to generalize expressions for g\ which can be used in the full region 
A © B © C © D. Let us first discuss the perturbative description of g\ in the region D. 
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The interpolation expression for g^ s valid at large x and small Q 2 is 



g™^, Q2) = (#2) y ^ £ (i) C^MSqM&ZZ*™ , (180) 

where a; = x + z, y = ln[(Q 2 + fj, 2 )//! 2 ] and the combined coefficient function C^L, and anomalous dimension h^JL b 
are given in Eq. (|161[) . 

Similarly, combining the shift and Eq. (|166[) . the expression for the singlet 171 in region D is: 

9?&*, Q 2 ) = IJZ^ (l)!^^ 6 ^"' 5 + <£U»)^-»)rt&>) + (181) 



19 qg 



Actually, Eqs. (|180ll8ip represent g\ not only in region D but also in the full region A © B © C © D. Indeed, they 
are expressed in terms of the shifted variables x, y and therefore can be used at any values of Q 2 . Also they include 
the total resummation of the leading logarithms, so can be used in the regions B and C. Finally, they contain non- 
logarithmic one-loop contributions 6 to the coefficient functions and anomalous dimensions obtained by the DGLAP- 
expressions, so this makes it possible to use them in the large- a; regions A and D. Let us remind that the use of 
the shift in Q 2 drives us out of the logarithmic accuracy and also recall our suggestion is to use different values of fi 
for the singlet and non-singlet components of <?i, given in Eqs. (|130I137[) . We proceed now to the applications of the 
results on g\ obtained so far.. 



XVII. PERTURB ATIVE Q 2 -POWER CORRECTIONS 

In this Section we discuss the power-l/(Q 2 ) fc corrections to g\ . Basically, there are various sources of such corrections 
but we focus only on those arising when the experimental results of g^ xp are confronted to the theoretical predictions 
g{ heor . The numerical analysis of the discrepancy between the non-singlet component of g^ xp and g t 1 heor shows (see 
for example Ref. [45| and Refs. therein) that 



( 9 n exp - (*rr or ~ fe E ^ m 

and the reason of the discrepancy is attributed to the impact of higher twists. Conventionally, the DGLAP expression 
of Eq. (|2U)) is used for describing (gi S ) & ° T and this is called the leading twist contribution. Naively one could expect 
from Eq. (|182p that the impact of the power corrections should increase when Q 2 decreases, especially in the limit 
Q 2 — > 0. On the contrary, the power corrections become negligible when Q 2 decreases down to values ~ 1 GeV 2 . 
We are going now to explain this behavior and give an alternative description of the power corrections. To this 
aim, first let us notice that the kinematic region of g^ s studied in Ref. [4|| mainly coincides with the region B © C 
and therefore the total resummation of the leading logarithms together with the shift of Q 2 should be included into 
expressions for g{ heor . Eq. (|179[) contains both terms and therefore in the region B © C 

where w = 2pq and /lijvs defined in Eq. (|130|) . The terms with Q 2 + [i 2 N s in Eq. (|183p can be expanded in the region 
B, where by definition Q 2 > n 2 NS , as follows: 



k=l 



(184) 



The second-loop contributions can be included similarly. 
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with 



^ns _ {-u + h NS )(-uJ + h NS - l)..(-w + /ijvs 



fc + 1) 



(185) 



Obviously, the power terms in the series of Eq. (] 1 84[) have a perturbative origin and have nothing to do with the higher 
twists. Such terms are absent in the Standard Approach. Of course we are aware that higher twists can contribute to 
gi but we argue that the perturbative power contributions of Eq. (|184[1 should be accounted for first, and only after 
a reliable estimate of the impact of the higher twists can be made. In contrast, in the region C where Q 2 < n%g, the 
power Q 2 -expansion takes the different form: 



^Q 2 



Q 2 



Mats 



G)"hEr«.>(£) 



*:=i 



(186) 



The power series in Eqs. (|184I186|) for large and small Q 2 are derived from the same formulae. However after the 
expansion has been made, they cannot be related to each other by simply varying Q 2 . Our estimate /ins ~ 1 GeV 
gives a natural explanation to the observation made in Refs. [45| that the power Q 2 -corrections die out when Q 2 
approaches values ~ 1 GeV 2 and do not appear at smaller values of Q 2 . Let us remind that our estimate of uns m 
Eq. (|130|) was obtained by studying the asymptotic properties of g^ s , i.e. absolutely independently of any analysis 
of the power corrections. We suggest that the new source of the power contributions given by Eqs. (|184I186[) can 
sizably change the conventional analysis of the higher twists contributions to the Polarized DIS. Obviously, the power 
expansion of the singlet g\ can be made quite similarly. 



XVIII. APPLICATION TO THE COMPASS EXPERIMENT 



Now let us discuss the application of our results to the recent COMPASS data on the singlet g\ . We consider here 
the results of Refs. fl6l |48||. The COMPASS experiment, carried out at the Super Proton Synchrotron at CERN. 
has investigated g± by measuring the asymmetries in the scattering of a polarized 160 GeV fi + -beam on polarized 
deuterons from a fixed 6 LiD target (see Ref. [Hj]). As there is only one target, the COMPASS collaboration can 
measure the singlet g\ only. Values of Q 2 at the COMPASS data are basically small: events with Q 2 < 1 GeV 2 
correspond to about 90% of the total data set. From Refs. one can conclude that the COMPASS kinematic 

region for measuring g x , Gcompass, is 



G 



COMPASS 



10 < x < 10" 



10" 1 GeV 2 < Q 2 < 1 GeV 2 



(187) 



This makes clear that the Standard Approach cannot be used for the analysis of the COMPASS data. On the contrary, 
our expressions (|1161 1181]) can be used in the COMPASS kinematic region. In the region Gcompass Q 2 *C /i| and 
therefore in this region 



-Q 2 



/4 



Q 2 



/'I 



<±) 



i 



k=l 



with 



( ±) _ j-uj + n±)(-uj + o± - i)..(-w + n± - k - 



k\ 

Substituting Eq. (|188[) into Eq. p8ip leads to the following expression: 

dgi(z,x,Q 2 ) 



gi {x,z,Q z ) ^gi(z) + (Q 2 /tii)- 



+ o((g 2 /4) 2 ^ 



dQ 2 /n 2 s 



(188) 



(189) 



(190) 



where z = fJ^/w. The first term in Eq. (|190p is 



<ko / 1 



2m V z 



C q {uj)8q + C g {uj)5g 



(191) 



As stated earlier, the combined coefficient functions C q%g include the total resummation of the leading logarithms of 
z and the non-logarithmic contributions ~ a s . They are defined as follows: 



G a — C„ + C, 



DGLAP 



Cq — C„ + C i 



DGLAP 



AC 



(192) 



3G 



The terms C° GLAP and C° GLAP in Eq. ([192]) are the NLO DGLAP coefficient functions and 



u{u - H gg ) -I^M 

H H - H H ' 9 ~~ 



° 9 c 2 _ u{Hgg + Hqq ) + H qq H gg H qg H gq ' 1 + W 2 ' 



£i . ~~y<j AC ^ 

'' " w 2 - ^(H gs + H qq ) + ff 9? ff gs - ff gg ff 39 ' 9 ~ lu 2 



The presence of the terms C®£ LAP in Eqs. (|192I193|) . as the DGLAP coefficient functions, may sound irrelevant 
or strange because the DGLAP description of g\ cannot be used in the small-Q 2 kinematics we are discussing. 
Nevertheless, the direct calculation of gi(z) to order ~ a s yields a contributions coinciding with the NLO DGLAP 
coefficient functions. In other words the presence of the coefficient functions has nothing to do with the Q 2 -evolution. 
Eq. (|190p explicitly shows that the Q 2 - dependence of gf in the region Gcompass should be weak, and also that 
gf practically does not depend on x in the kinematical region Gcompass , even at very small x. However, its 
absolute value cannot be fixed from theoretical grounds. Indeed, Eq. (|191[) implies that sign of gf (z) at any given z 
is determined by the interplay between the quark and gluon contributions and eventually depends on Sq/Sg, which 
cannot be determined theoretically. At the same time, Eq. (|190l) predicts that the z-dependence of gf is pretty 
far from being trivial, so the experimental investigation of this dependence would be quite interesting: it can yield 
information about the initial quark and gluon densities. The predictions of the essential independence of gf on x 
made in Ref. 16], was confirmed in Ref. where a flat dependence of gf was found. More precisely, Ref. 
reported that 

gf » (194) 

in the region Gcompass with small errors. Unfortunately, the COMPASS data do not allow to study the z-dependence 
of gf in the proper way. Nevertheless, the COMPASS result (I194j) was used in Ref. [48| in order to obtain some rough 
estimates for dq/Sg. Below we consider this issue in detail. 

A. Interpretation of the COMPASS data on gf 

The variable w — 2pq in the COMPASS experiment runs in the interval 

30 GeV 2 <w< 270 GeV 2 (195) 

and therefore the range for z is 

l<z<0.1. (196) 
The variable z is related to the standard variable v = w/(2M) measured in GeV, with M = 1 GeV: 

\2MJv v 1 v ' 

so the region (|196[) covered in the COMPASS experiment corresponds to the i/-region (in GeV) 

15<^< 150. (198) 

We remind that only the x -dependence of gf was studied in the COMPASS experiment. The values of w and Q 2 
were not reported in the COMPASS data, which makes impossible the straightforward application of Eq. (|191[) to the 
COMPASS results. However there are several options for the interpretation of Eq. (|194[) and below we consider them 
in detail: 
Option (i): 

Eq. (|194[) means that gf(z) = for any z from the whole interval of Eq. (|196p . 
In this case Eq. (|191[) implies a strong correlation between 5q and 5g at any lu : 

C q (uj)Sq(uj) + C g {Lu)5g(cu) = 0. (199) 



We don't find theoretical grounds for understanding this fact and think that next option is more realistic. 
Option (ii): 



37 



Eq. (|194p holds in the average, namely in the region (|196|) : 

< gf{z) >= . 



(200) 



Obviously, in order to fulfill the Eq. <|200[) . g\(z) should acquire both positive and negative values in the region (|196|) . 
This could be realized by an appropriate choice for the initial parton densities Sq(z) and 5g(z). In ref. [l4| we suggested 
that in region B one can approximate the initial parton densities by constants. Guided by this result, we suggested 
in Ref. [16| to approximate Sq(z) and Sg(z) at small z by simple constants to get a rough estimate. However, in the 
COMPASS region (|196p z is not small enough to use such a simple approximation. As the DGLAP-fits from Ref. [j| 
work quite well and also other parameterizations have a similar structure, we suggest a similar but regular fit : 



The main difference with Ref. [5j is in the absence of the power factors z a while the terms in the brackets in 
Eq. (|201|) and in Ref. [f| coincide ( x in Ref. || is replaced by z in Eq. (|20ip ). Indeed in Ref. as also discussed 
previously, we have proved that the role played by the singular terms x~ a in the DGLAP fits is to mimic the total 
resummation of ln fe (l/x) . On the other hand, we would like to keep the same ratio Sq/Sg as in Ref. Q and therefore 
we also change the power factor for Sq in Eq. (1201 j) . Now it is easy to check that the fit (|201|) do not lead to a flat 
z-dependence for 171 and cannot keep gi(z) — in the whole COMPASS region (|187j) . 

In more detail by substitution of Eq. (1201)) into Eq. (|181[) and performing the integration over uj numerically, with 
fixed and positive N q , and varying the values of N g , we plot our results in Fig. [11] . 

By a close inspection of the various configurations shown, we can easily conclude that these fits could be compatible 
with Eq. (p00| only if N g > and N g > N q . 

As the way of averaging g\ over z in the COMPASS data is unknown to us, we can try another possibility, 
approximating 



where < z >= 0.25 (i.e. < v >« 60 GeV) is the mean value of z from the region (|196p . Then using Eqs. ()181|201j) . 
keeping positive N q and varying N g , as shown in Fig. [TTJ we suggest again that N g are positive and N g > N q . 



To conclude this Sect., let us make a brief comment on the parametrization of 171. In the Born approximation gi 
is given by Eq. ((J) and depends on the only argument x which corresponds to the famous scaling in the DIS. The 
radiative corrections in the higher loops bring the violation of scaling, so g\ acquires, additionally, the Q 2 -dependence. 
At this stage one can parameterize g\ by the set of variables x, Q 2 or, alternatively, w, Q 2 , or v, Q 2 . As it is well- 
known, in the target rest frame w = 2pq — 2M(E — E') where M is the target (nucleon) mass and E (£") is the 
energy of the incoming (outgoing) lepton. Then Q 2 in the same frame involves the above energies and the scattering 
angle 8. Therefore, both x and Q 2 depend on E' and 9 and these variables are not always independent of each other. 
Indeed, there are experiments where the x -dependence of g± is measured at fixed Q 2 ; then Q 2 is varied to another 
value and the x -dependence is studied again. In this case x and Q 2 are really independent variables. In the opposite 
case of fixed w and varying Q 2 , Then x = x(Q 2 ) and therefore these variables are not exactly independent, so it is 
more convenient to use w instead of x, as independent variable. These examples show that using w and Q 2 instead 
of the standard set x, Q 2 could be more convenient for gi . In particular, when Q 2 is very small, using x instead 
of w becomes really inconvenient. Nevertheless, the w -dependence of g\ in the COMPASS experiments predicted in 
Eq. (|190p can be clearly extracted from the dependence on x and Q 2 . Indeed, the variable x defined in Eq. (|178[) can 
be written in the following form: 



8q{z) = N q z{l ~ z) 3 (l + 3z), Sg(z) = N g (l - z) 4 (l + 3z) . 



(201) 



<gx(z) >« gi (<z>) = , 



(202) 



B. Comments on the measurement of gi in different kinematic regions. 



x = /i 2 x/Q 2 + x fa \j?xjQ 2 . 



(203) 



Although gi depends on x and on Q 2 at small Q 2 very weakly, its dependence on x/Q 2 is quite essential. Indeed 
Eq. (|190p can be regarded as a sort of a new scaling law where g\ depends on the variable x/Q 2 only. 
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FIG. 11: The v -dependence of gi(y), with 8q, Sg defined in Eq. (|201[) , for N q — 0.5 and different values of N g : (a) -1.5, (b) 
-0.5, (c) 0, (d) 0.5, (e) 2, (g) 3.5; the COMPASS v -region corresponds to Eq. JT98J . 
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XIX. SUMMARY AND OUTLOOK 

In the present paper we have presented an overview of our results on the spin structure function g\ at arbitrary 
x and Q 2 . We have divided the whole kinematic region of x and Q 2 into the set of four regions A - D defined in 
Eqs. (|51 ITDTfT5|) and considered gi in each of these regions. The region A is covered by the Standard Approach, based 
on the DGLAP evolution equations. This is briefly discussed in Sect. II. The application of the integral transforms to 
gx is given in Sect. III. In Sect. IV we have discussed in detail the parametrization of a s and shown that the popular 
parametrization a s — a s (k 2 _) is valid at large x only. Otherwise it should be replaced by the effective coupling a e Jf 
defined in Eq. (|^|) . When /z 2 obeys Eq. (|5D|) . a e Jf can be approximated by a s {k 2 L / (3) (with (3 being the longitudinal 
Sudakov variable) and when, in addition, x is large, it can be simplified down to a s (k 2 _). According to Eq. (|50p . the 
deviation of ofj' from a s strongly depends on /i 2 . For example, when /i 2 = 2.5 GeV 2 , u s (n 2 ) / 'a e J '•' (/J 2 ) ~ 0.9 but 
very quickly a s (ix 2 )/a e J ' (/i 2 ) ~ 0.5 at /i 2 sw 1 GeV 2 , which is a typical DGLAP starting point of the Q 2 -evolution 
The small- a; region B, where the total resummation of the leading logarithms of x is essential, was considered in 
Sect. V. We account there for the resummation of the leading logarithms by solving an Infrared Evolution Equations 
(IREE), so in Sect. V the essence of the method together with the technology of IREE was discussed and the IREE 
for gi were obtained. These IREE involve new anomalous dimensions and coefficient functions. Explicit expressions 
for them were obtained in Sect. VI. Then the expression of Eq. (| L06[) for the non-singlet component of g\ in the region 
B was obtained in Sect. VII, while in Sect. VIII the result for the singlet g\ in region B was given in Eq. (| L L6|) . 
Obviously, the impact of the total resummation of logarithms of x is big at small x and becomes maximal at x — > 0, 
where the expressions in Eqs. (|106|116|) behave asymptotically as in Eqs. (|124|134j) . The small-x asymptotics of 
the non-singlet g\ is considered in detail in Sect. IX while Sect. X contains the asymptotics of the singlet g\. Both 
asymptotic behaviors are of the Regge type. The estimates for their the intercepts are given by Eqs. (| L 29|) and (]L36|) . 
The small-x rise of g\ predicted by Eqs. (|124I134[) is much steeper than the well-known small- a; DGLAP prediction 
in Eq. (|132[) . On the other hand, a numerical analysis shows that the use of the asymptotic formulae at the presently 
available x is not reliable, so Eqs. (|106I116[) should not be replaced by their asymptotic expressions of Eqs. (|124I134[) 
in the region B. The comparison of our results to the DGLAP -expressions for gi, which is impossible without fixing 
the initial parton densities Sq and 5g, shows that the impact of the resummation of the logarithms becomes essential 
for values smaller x « 10 -2 In Sect. XII, by considering in detail a standard DGLAP fit (f24"f for the initial parton 
densities, we have shown that the singular factors in the fits mimic the total resummation of the logarithms and 
provide the rise of g\ at small x which is observed in the experimental data. On the other hand, when the total 
resummation of the logarithms is taken into account, the singular factors in the fits can be dropped. This allows one 
to simplify the parametrization of the initial parton densities.. 

The Reggeon structure in the two approaches has been discussed in Sect. XIII, However, in the case of the SA those 
Reggeons are, in a sense, fictitious: they are generated by the fits for the initial parton densities and because of that 
are present at any x instead of appearing in the asymptotic expressions at x — > 0. 

The total resummation of the small-x logarithms is important in the region B, but also non-logarithmic contributions 
are quite essential in the DGLAP region A, where they are accounted for to NLO accuracy. Then the interpolating 
expressions for g\ are presented in Eqs. (|163p and (|166[) . On one hand, they almost coincide with gDGLAP j n ^ e 
region A and on the other hand, with Eqs. (|106I116[) in the region B, and at the same time do not require the use of 
the singular parameterizations of the parton densities. 

The small-Q 2 region C is absolutely beyond the reach of DGLAP. On the other hand, the analysis of the Feynman 
diagrams contributing to gi shows that a shift of Q 2 allows us to extend Eqs. (|106I116[) into the region C. Similarly, 
Eqs. (|163ll66p can be extended into the region D. Eventually we arrived at the Eqs. (|180I181|) which are the interpo- 
lation expressions which can describe g\ in the whole region A © B © C © D. This was the subject of Sects. XIV- XVI. 
This shift of Q 2 , given in Eq. (|167p . inevitably causes the appearance of power 1/Q 2 -corrections. Q 2 power corrections 
were found earlier phenomenologically by confronting gDGLAP an( j experimental data. They were attributed to 
the impact of higher twists. In Sect. XVII we argued that the role played by the higher twists can be estimated 
reliably only after accounting for the pertubative power corrections. 

Finally, in Sect. XVIII we have used the small-Q 2 description of the singlet g\ in Eq. (|18ip for the interpretation of 
the recent COMPASS data. First we have shown that g\ in the COMPASS kinematic region does not depend on x, 
even at very small x. Then, we have suggested that the COMPASS data are compatible with positive gluon densities 
. We also argued in favor of studying the dependence of g\ on 2pq in the COMPASS experiment rather than on x, in 
order to estimate the ratio Sg/5q. 
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APPENDIX A: SIMPLIFICATION OF THE COLOR STRUCTURE OF THE FORWARD SCATTERING 

AMPLITUDES 

We consider below in more detail the color structure of the Born amplitude A Born defined in Eq. (I42[) . As the 
external partons in A Born are quarks or gluons, 

A Born = _ C (col) 4 "(-jfrb/^CPlV '(Plb/X |H*) ( M) 

s + ie 1 ' 

where the SU(3) -matrix C^ co ^ describes the color structure of A Born . When all external partons are quarks, 
(j( c ° l ) — t a t a , with t a (a — 1, .., 8) being the SU(3) -generators in the fundamental (three-dimensional) representation; 
when the quarks are replaced by gluons, t a are replaced by the SU(3) -generators T a in the vector representation. 
Each of the initial and final color two-parton states in Eq. (|A1[) corresponds to a reducible representation of SU (3) 
and can be expanded into a sum of irreducible states. It is convenient to do it in the t -channel where the amplitude 
A Born describes the quark- antiquark annihilation qq — > q'q' and therefore the irreducible initial qi<p and final q p q q 
color states (with i,j,p, q = 1, 2, 3) are the singlet (5) and octet (V). The initial 3 (g> 3 color state in the t -channel 
state is qiq 3 . It can be expanded into the sum of the singlet and octet, each one is the irreducible state: 3® 3 = 1 08. 
We denote them (qiq J )s and (qiq J )y respectively. It can be done by applying the projection operators P$ and Py to 
the quark-antiquark states: 

(m 3 )s = (Psifa'?' , (qi^)v = (iV)£W (A2) 

where 

p S = ^'4' = ^ a )i{t% ■ ( A3 ) 

Obviously, these operators are orthogonal to each other and the factors 1/N (with N = 3) and 2 in Eq. (|A3|) are intro- 
duced to guarantee the property P 2 = P for each of Ps and Py ■ Let us notice that | \P$ \ | 2 = Tr[PgPs] — 1, | \Pv 1 1 2 = 
Tr[Py Py] = 2NCf where Cp = (N 2 — 1)/2N. Applying the projection operators Ps and Py to the color factor 
C(coi) in Eq . allows 

us to write down C^ co1 ^ as the sum of the scalar Cg CO '' ) and octet Cy' color factors: 

(c ic ol)) rp = C ^ l \p s )f q + C^ l) (Py)% , (A4) 

where 

r {col) _ ( P g)jp( ta )i( ta )g _ 1 Trh a, a] _ r 

5 " W -jjTr[tt\-C Fl 



v - \\Py\\ 2 ~ 2NC F ~ 2N ■ ( ' 



lM ) _ (Pv)%(t a ) j(r) p q _ 2Tr[tnHH b ) _ 1 
' V ~ \\Pv\\ 
Substituting Eq. (| A4|) into Eq. (|A1|) . we rewrite it as 

A Bom = A Bomp s + A Born Py , (A6) 

with 

A Born _ _(7( co/ ) 47rQ , u(-Plh» U (Pl) U '(Plhli U '(-P2) 

s s s + ie 
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When the external quarks in Eq. (|A1|) are replaced by gluons, generators t a in the factor C^ ^ are replaced by the 
SU(3) -generators T a in the vector representation. It allows one to generalize Eq. (|A7[) to the gluons case, expanding 
the initial t -channel gluon state 8®8 into scalar and octet (see Refs. fiol ITTj ] for detail). The projection operators 
Pg 99 ^ and Py 9%> projecting the two-gluon t -channel state on the scalar and octet states are: 

P ( s 99) = J^I—J <W<U , Pl? 9) = ^(T c ) a , b/ (T c ) ab . (A8) 

Strictly speaking, the normalization for projector operators with the gluon-quark transitions can be arbitrary but in 
order to match DGLAP it can be chosen as follows: 



x(t c )i(T c ) ab , (A9) 

In Eqs. (|A8|A9j) we have kept the notation i,j for the quark color states while a, b, a', b' denote the gluon states. In 
contrast to the quark-quark case, the expansion of 8(g) 8 into the sum of the irreducible SU(3) -representations includes 
the singlet 1, the antisymmetric 8a and symmetric 8s octets and other contributions which cannot be organized out 
of the gluon fields and therefore can be left out, so for the gluons 8®8 = lffi8 A ffi8 s . The symmetric octet 8 S does 
not contribute to gi with the leading logarithmic accuracy (see Ref. [lOl ] for detail). An additional argument in favor 
of neglecting the amplitudes with high color dimensions is that they die out quickly with energy. 



p(<?s) _ 



1 



p, 



(91) 



N 2 
1 

AT 



5U ab 



p(<79) _ 



3(99) 



APPENDIX B: NON-SINGLET CONTRIBUTION TO THE STRUCTURE FUNCTION F x 

The technology for calculating the leading logarithmic contributions to F^ s (x, Q 2 ) is quite similar to the one for 
g^ s (x,Q 2 ) (see Refs. [ll| for detail). Similarly to gi S , F{^ s is expressed through the forward Compton amplitude 
T^g in the following way: 

F» s = ±ST£j. (Bl) 

The superscript (+) in Eq. (|B 1 1) stands for the positive signature. Then it is convenient to define the Mellin amplitude 
Ff s (uj, y) related to through the Mellin transform ([53]). As a consequence, the IREE for Ff s is almost identical 
to Eq. !p]): 



U F» S { U , V ) + 1 d ^ y > = (i + V) £ff(o;)*n«,v) ■ (B2) 

The amplitude Lgg (ui) again corresponds to the quark-quark scattering but its signature is now positive. It should be 
found independently. Eq. (|B2p can be solved similarly to Eq. The only difference between them is replacement 

of tiNS by h( + \ The amplitude obeys the following IREE: 



^(+»=t(+) + (l + V)(hW) 2 . (B3) 

The difference between Eq. (|B3[) and Eq. ([80]) is the inhomogeneous term M + ) . It is expressed through a qq defined in 
Eq. and V q ( + ] : 

b^ = a qq + Vj+\ (B4) 

Similarly to V qq introduced in Eq. (|88p . V qq ^ is expressed in terms of m qq defined in Eq. (|89[) and a new quantity 
instead of D: 



99 



m qq D^ 



(B5) 
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with 



D(u) 



2b 2 



dpe-^f In ((p + r])/ri) 



p + i] 



1 



(p + rj) 2 + n 2 p + i] 



So, the expression for F{ b in region B is 



(B6) 



where 



2 />200 



du) 



2 ./_„„„ 27TI 



C 



(+) 

NS 



2uj 



(B7) 



(B8) 



and 



ft<+>(u;) = (1/2) L - y/u2-B$(w) 



(B9) 



B<+J = 4(1 + A M w)&&> . (BIO) 



The small- x asymptotics of F( fS (x,Q 2 ) is also of the Regge type but the value of the intercept A^J is smaller 
than the one of gf 18 : 



A ( + ] s = 0.38 



(Bll) 



APPENDIX C: CONVOLUTION OF TWO AMPLITUDES 



(Pi) 



Let us consider the t- channel convolution of two amplitudes: Q = Af 1 Aif 2 of the scattering amplitudes A ( { 
and A^ 2 ^ where p\.2 = ± stand for the signatures. It is convenient to describe amplitudes M^"^ in terms of the 
invariant amplitudes M a . For example, the invariant amplitude for the quark-antiquark forward annihilation 

l{pi) + 9(^2) — * q(p'i) + are introduced as follows: 



A {±) = Mil M (±) 
s 



(CI) 



where j„ are the quark currents. We remind we use the Feynman gauge. 

It is also convenient to use the asymptotics of the Sommerfeld- Watson transform (often called the Mellin represen- 
tation) for each of those amplitudes in the following form: 



with a = 1,2 . The signature factor 
and the transform inverse to Eq. (|C2I) is 



±1 (1 ± 1) + 



ttuj J 2 s \p^y 



1 f 00 ds 1 s \ %M ± 3 U M 
m 2 



(C2) 



(C3) 



(C4) 



Using Eqs. (|C1IC2[) and skipping the overall factor j v j v / ' s allows us to write the convolution Q through the invariant 
convolution Qi nv as 



Q 



Juju , 



)PlP2 

' inv 



(C5) 
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with 7 

where the factor 2k\ appears as the result of simplifying the spinor structure, si = 2p\k, S2 — 2p2k. Both of them 
are understood as Si t 2 + tt. 

For integration over k in Eq. (|C6|) we use the Sudakov variables (fT5|) : 

« 2 = I^I ^^^ ) (^)^ ) (^)/i W (^)/ a W (^)*^* ) (^.^.-/M a ) (07) 

with 

•*"■>(«.«../*■> = / ^ 2 ^ ^-| + .^ (m)"'(m)"^A ■ (C8) 

Let us first integrate Eq. (|C8j) over a. The integration can be done in the complex plane by applying the Cauchy 
formula. The singularities in the complex a -plane are the double pole sa/3 — k\ + ie = and the cut from the Mellin 
factor (sa) Wl . The integration yields a non-zero result when the pole and the cut have opposite imaginary parts. 
The imaginary part of the cut is positive while the imaginary part of the pole is negative provided f3 > 0. When 
/3 < 0, both singularities have positive imaginary parts and therefore the integration over a yields zero. Closing up 
the integration contour in the lower hemi-plane and taking the residue of the pole 

a = (k 2 ± -ie)/s/3 (C9) 

we perform the integration over a and arrive at 




(CIO) 



Therefore we obtain the following expression for Qi nv : 

«r = A r pp****** ******** [±(j,y -±(±f] . (cu , 

sir'' J_ loc 2m 2m u>2 — wi lu>2 V/x / uii v/z / J 

Eq. (|C11[) involves two integrations, however one of them can be done easily. The integration lines over in 
Eq. (|C11[) are parallel to the imaginary axes and lie to the right of the rightmost singularities of f\$. Let us assume 
that additionally to it 



< 5Rwi < nuj 2 . (C12) 

The opposite case Sftwi > $iu>2 can be discussed similarly. The integrand of Eq. (jClip includes two similar terms in 
the squared brackets. Let us focus on integrating the first term, (s//! 2 ) 2 and let us integrate this part of Eq. (IClip 
with respect to u)\. In this case closing up the u)\ -integration contour to the left involves accounting for singularities 
of fi(uj\). On the contrary, when we close up the lo\ -contour to the right, the only singularity inside the contour is 
the pole 1/(cl>2 — uii), so we can do this integration without considering f\, just by taking the residue at u>i — u>2- At 
the same time, integrating the remaining, proportional to (s/fj, 2 ) 1 part of Eq. (|C11[) with respect to uj\ yields zero. 
Therefore Eq. (|C11[) is reduced to the simpler form: 



QKT - r ^(VmT^M^V) rt Pl){uj) f 2){Lj) (ci3) 



7 The factor % is the product of the overall factor — % and (±i) 2 from the t-channel quark or gluon propagators. 
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Obviously, 

^>M^ 2 )H » (1/4) [(1 + PiP 2 ) - zvr W (P 1 + P 2 )] = <W 2 (C14) 

where Pi 2 — ±1 ■ It means that the leading contribution to Qf,^ 2 is diagonal in the signatures. Strictly speaking, 
this should be checked in advance, with using the Sommcrfeld- Watson transform where analytical properties of the 
involved amplitudes are explicitly accounted for. 



APPENDIX D: OFF-SHELL INVARIANT AMPLITUDE M IN EQ. ([169)) 



The invariant amplitude M in Eq. (|169|) is a generic notation for the following invariant off-shell am- 
plitudes: M NS (2pk,k 2 ) contributing to the Bethe-Salpeter for g^ s and the flavor singlet amplitudes 
M qq (2pk, k 2 ), M qg (2pk, k 2 ). They are related to the amplitudes A qq (2pk,k 2 ), A qg (2pk,k 2 ) of the forward quark- 
quark and quark-gluon scattering: 



A„ 



u(-k)j\u(p)u(p)j\u(-k) 



(p + k) 2 



ic 



M qq {2pk,k 2 ) , 



(Dl) 



, u{ -k)-fx(p - k)j u u(-k) 
(p + k) 2 + it 



M qg (2pk,k 2 ) 



The outgoing momenta —A; in Eq. (|D1|) are assigned to the final (upper) off-shell quarks whereas the initial quarks 
and gluons have momenta p ; e\ and are the polarization vectors of the gluons. We remind that the amplitudes 
M qq (2pk,k 2 ), M qg (2pk,k 2 ) are off-shell and therefore they differ from the amplitudes M ik introduced in Eq. (f73"j) : 
M qq and M qg logarithmically depend on two arguments, 2pk and k 2 , i.e. 



M, 



gq 



M qq (p, z); M qg = M qj g(p, z) 



where p = \n{2pk/ p 2 ), z — \ii{k 2 / p 2 ) . 

It is convenient to introduce the Mellin amplitudes tpik{uj,z) conjugate to Mik(p,z) through Eq. 
for ipik(oJ, z) is quite similar to Eqs. (|75l76p . and we put here all = for the sake of simplicity: 



(D2) 
The IREE 



dtp 



NS 



dz 
dz 

dt Pqg 



LOtp qq 
-UJip qg 



NS 



1 

" 8^ 

1 

In- 
1 



V NS L qq (u;) , 



^i^qqLqqiu) + ^<PqgL gq (uj) , 
1 



(D3) 



The amplitudes L NS , in Eq. (|D3[) are on-shell, so in accordance with Eq. ([77)) they can be expressed in terms 
of h NS (uj) = (1/8tt 2 )L ns (ui), h ik {uj) = (l/8vr 2 )L ifc (w) which are obtained in Eqs. (|104I95|) . General solutions to the 
linear equations (|E>3|) can easily be found (cf Eq. (|107p ): 



ip qq = * 1 (w)e- wz+z °(+) + * 2 Me" a "• ' 2: 



(D4) 



tp qg = *i(w) 



X 



R 



2h n 



X-y/R 

2hnn 



with f2(±), X and R defined in Eqs. (|110H108I109[) respectively whereas 'J'i^ should be specified. Therefore, E(2pk.k 2 ) 
used in Eq. (I168j) can be any of E NS {2pk.k 2 ) 1 E qq (2pk.k 2 ), E qg (2pk.k 2 ) given by the following expressions: 



E NS (2pk.k 2 ) = 
E qq [2pk.k 2 ) = 
E qg {2pk.k 2 ) = 



^Lf^Y^S^z,^ 



—zoo 



2m V k 2 J 



(D5) 
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duj /2pk\ L 
2m\W) 

2n~i\~k~ 2 ~) 



*i(w)e z °<+> +* 2 (^)e z °(-' 

x + Vr 



2h„ 



2h n 
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In order to specify ^i^, we use the obvious matching condition: 

(p NS (ui, z = 0) = 8n 2 h NS (uj) , f qq {uJ, z = 0) = 8Tr 2 h qq (u) , V qg (u, = 0) = 8ir 2 h qg {uj) . (D6) 
It immediately fixes ip NS : 

<p"S = 87r 2^S (w)e ,[- W +^ S (c)] (D7) 

and leads to the explicit expressions for ^1,2: 

^ = 0,2 [ 2h qg h gq - h gq( h gg ~ fe w ~ V^)] ^ = 8?T 2 [ ~ 2h qg h gq + fe w(V? ~ ^gg + V^)] / Dg N 

1 77 2v^R ' 2 2v^R 

APPENDIX E: CALCULATING THE SMALL-i ASYMPTOTICS OF THE NON-SINGLET g x . 

Eq. (| 1 06[) for s l(a;, Q 2 ) in region B can be written as follows: 

9 r^Q 2 )= r (ei) 

J-lOO 27TI 

where the phase <J> is 

2/ 



$(w,x,Q 2 ) =wln(l/:r)+lii(7 i vsH + yft J vsH =w£ + hiC N s(u) - ^lu 2 - B ns {lu) (E2) 



where we have used the expression (| 104|) and denoted £ = ln(w/ y/ Q 2 p 2 ). We remind that y — ln(Q 2 /p 2 ). We are 
going to calculate the asymptotics of g^ s at x — * and fixed Q 2 , i.e. at fixed Q 2 and w — ► 00. The standard way to 
calculate asymptotics is to apply the saddle-point method to Eq. (|E1[) . According to it, 

g? S ~ H^^o, Q2 )e *o(-o,»,Q 2 ) j (E3) 

with the stationary phase = ^K^o, w j Q 2 ) and the stationary point ujq is defined from the requirement d^/duj = 0, 
i.e. loq is a solution to the equation 

CWsM 4 ^2 _ BjVS ( w ) 

Obviously, Eq. (|E4|) can have much more the one solution. In this case luq is the solution with the largest SRw. 
It is often called the rightmost stationary point. Substituting the explicit expressions (|105[) for Cns, we transform 
Eq. JEH into 



faB NS \/u 2 -B NS = (Bns - cjB' ns /2)(uj - y/u; 2 - B NS ) + ^uB NS {u - B' NS /2) . (E5) 

Obviously this equation cannot be solved analytically. The analytical solution can be found for the particular case 
when B does not depend on u> (it corresponds to the case of fixed a s ) and y = 0. In this case Eq. (|E4|) can be reduced 
to the algebraic equation 

+ (2/p)oj 3 ~ uo 2 B NS - {2B NS /p)w - (B NS /p 2 ) = (E6) 

where p = ln(w/p 2 ). Eq. (|E6|) has four roots which can be found with using the known from the literature Ferrari 
formulae but only two of them, namely lu = ±^JBns do not go to zero when w — > 00. Obviously, in this case the 
rightmost root is uio = \/Bms- It is easy to make this conclusion, without solving Eq. (|E6[) . Indeed, Eq. (|E6|) can be 
written as 

u 2 {u 2 - B NS ) + {2/ P )lu(lo 2 - B NS ) - (B NS /p 2 ) = . (E7) 

When the terms ~ 1/p 2 and ~ 1/p are dropped, Eq. (|E7j) can be solved immediately and the rightmost root can 
easily be found. Applying the same arguments allows one to solve Eq. (|E4|) at w — > drives us to conclude that the 
rightmost and non-vanishing at w — > 00 root of Eq. (|E4|) does not depend on y and it can be found as the rightmost 
root of the much simpler equation uj 2 = Bjys(u) as is stated in Eq. (|125|) . This leads to the Regge asymptotics of 
Eq. (|124p quite different from the well-known DGLAP asymptotics ()132[) . Let us notice that, in the perfect agreement 
with the concepts of the phenomenological Regge theory, this root corresponds to the branching point singularity of 
Eq. jE4|. 
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APPENDIX F: THE DGLAP SMALL-i ASYMPTOTICS 



Let us remind how the DGLAP asymptotics of Eq. (|132[) was obtained. As this topic is well-known, for the 

sake of simplicity we consider 9i£>glap with the LO accuracy. When the singular term in Sq is absent, the small- a; 
asymptotics of 9idclap can a ^ so ^ e obtained with the saddle-point method. Similarly to Eq. (|E1[) . 9idql A p can 
be written as 



where the phase Qdglap is 



*dglap=wHV*) + I ^^4z^7 (0) , (F2) 



with y(°' being given by Eq. (f2"T]) . The bulk of the integral in Eq. (|F1[) comes from the region of to obeying 

wln(l/x)<l (F3) 

because the factor e"' 11 ' 1 / 1 ' in Eq. (|F1[) strongly oscillates beyond this region. So, the values of to mainly contributing 
to the integral become small when x — > 0. As a consequence, the most important term in is now the singular in 
to term Adglap(Q 2 )/u, with 



Adglap(Q 2 ) = / % 



dk\ a s (k 2 ± )C F 



2n 



(F4) 



Therefore, approximately 



®dglap ~ lo ]n(l/x) + A dglap {Q 2 )/lu . (F5) 
The equation for the stationary point is 

®'dglap - In(V^) - A dglap (Q 2 )/lu 2 = , (F6) 

which leads to the stationary point 



u° GLAP = ^A DGLAP (Q2)/ \n(l/x) (F7) 

and eventually to the DGLAP -asymptotic given in Eq. (|132|) . Contrary to the case considered in Appendix E, the 
DGLAP stationary point depends on ln(l/x). Eq. (|F7j) shows that the Q 2 -dependence in the DGLAP asymptotics 
follows from the DGLAP parametrization a s = a s (fc^) and mostly from keeping Q 2 as the upper limit of the 
integration in Eq. (|F4|) . The latter takes place because of the use the DGLAP ordering (|13p . However in the small- a; 
region the ordering ()13|) becomes unreliable and should be replaced by the ordering of Eq. (|14|) where the upper limit 
is w. Obviously, the Q 2 -dependence in Eq. (|F7[) vanishes after replacing Q 2 by w. We remind that the DGLAP 
asymptotics (1132[) can be obtained only under the assumption that the initial parton densities are not singular at 
x — > otherwise the asymptotics (|F7[) is changed for the Regge asymptotics (|149|) . 



APPENDIX G: THE SMALL-i ASYMPTOTICS OF g? b WITH THE TRUNCATED SERIES FOR THE 
COEFFICIENT FUNCTIONS AND ANOMALOUS DIMENSIONS. 



We consider here the case where the coefficient functions and anomalous dimensions are calculated in high orders 
in a s , however without the total resummation of those contributions. According to Eq. (|F3[) the essential values of to 
in Eq. (IF1|) are small at x — > 0, so the most important terms in expressions for the non-singlet coefficient functions 
and anomalous dimensions are the most singular terms in u>, i.e. the double-logarithmic contributions. They can 
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be obtained, expanding Eqs. (|104i I105|) into series absolutely in the same way as was done in Eqs. (|144i I145|) . The 
expansion for the coefficient function in expressions (|106[) and $2U§ have the same form 

however, with different a. For Cns 01 Eq. p06[) 

a = S ws H, (G2) 

with -Bats given by Eq. ([94]) . Alternatively, when the DL contributions to Cns are calculated in the DGLAP 
framework, a — cldglap'- 

cldglap = a s (Q 2 )C F /(2n) . (G3) 

In contrast, the expansions for the exponent in Eq. (fT7|) involves integrations of a s (k 2 _) and therefore those two cases 
look quite different. As this difference it is not essential for the topic we consider here, we will use the approximation 
of fixed QCD coupling for the DGLAP description of g^ s . In this case the most singular, i.e. DL contributions to the 
DGLAP expression for the anomalous dimension of g^ can easily be obtained from the series for Hns- By doing 
so, we arrive at the following series for the exponent in Eq. (|106[) : 



Y DL = yH NS = y 



n 

a a 



Auj 16w 3 32w 5 



(G4) 



with a defined in Eq. (|G2|) and y = \n(Q 2 / fj, 2 ) , whereas in the DGLAP case the DL contribution to the exponent in 
Eq. (fT7|) can again be obtained with replacement a by cldglap'- 



T DL I dk l~(,, n \ a DGLAPy a DGLApy 2 , _ a DGLApV^ , ,(,-, 
^ DGLAP = / _ "7 c 2 _ 7l w : a s) ^ r C 2 — h C 3 — h ... (Ut>J 



with C2, C3 being numerical factors. Therefore, in the n-th order in a s the most singular contribution to Cns can be 
written as follows: 

41 = C(„)^ + 0(1/^ 2 "- 1 ) (G6) 

where c/ n ) is a numerical factor. It can be obtained with further expansion of Eq. (| 1 05[) into series. Similarly, the 
most singular contributions, Y DL " and T^glap are 



r (n) DL _ - a " p(n) DL _ apGLApV" 1 (C7] 
~y n ,,2n-l ' DGLAP n . .« J ^VU 



with c„, c„ be numerical factors. 

The phase $ in the Mellin integrals (|EIj) and (jFljl can now be written as follows: 

a n 

$ « w( + ]nC NS +yc n ^— T , (G8) 

rft ~ , , Wl _i_ lr, r<DGLAP , „ a DGLApV n 

Qdglap ~ wln(l/a;) + InG^j + c„ , 

where we have denoted £ = (1/2) ln(Q 2 /(a; 2 /z 2 )). Let us first consider the asymptotics of g^ s at Q 2 ~ M 2 - I n this 
case the last term in each of the equations in (|G8[) is zero and the stationary point is determined from the following 
equation: 

C+^=0. (G9) 
Cns 

Obviously, this equation does not have solutions leading to the Regge behavior of g^ s because all terms in Eq. (|G1[) 
are positive. 
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Let us consider now the case of large Q 
are: 



2. 



) 2 3> M 2 - The equations for the stationary points of the phases in Eq. (|G8|) 



DGLAP 



= ln(l/x) + 



c (n) (2r 



yc 



n' DGLAP 
°jVS 

r'DGLAP 

ViVS 



,2n 



o , 



(GIO) 



a DGLApV _ q 



The use of Eq. (|G6p for the coefficient function makes easy to see that the term C' ns /Cns and C N g GLAP /C® GLAP 
are proportional to l/o>, so these terms are much less singular than the last terms in Eq. (|G10j) and therefore they 
can be dropped. After that solving Eq. (|G10|) is easy and we arrive at the following approximate expression for the 
stationary point: 



LOO 



((2n - l)c (n V n VC) 



l/2n 



.DGLAP 



(ncna n DGLAP y n I ln{l/x))) 



l/(n+l) 



(Gil) 



and leads to the following asymptotics: 



9? S 



exp 



(.<■>) 



l/2n 



IS 



l/2n Al-l/2n) 



NS DGLAP 

9i 



exp 



l/(n+l) 



l/(n+l) 



(ln(l/V 



(i-i/(»+i)) 



(G12) 

Obviously, Eq. (|G12|) coincides with the LO DGLAP asymptotics (jl32[) at n = 1. Eq. (|G12|1 demonstrates explicitly 
that the asymptotics of s DGLAP always depends on Q 2 and the Regge behavior of g^ s and g^ s DGLAP cannot 
be achieved at any fixed n. On the other hand, the Regge behavior of g^ s and g^ s DGLAP is approached closer and 
closer when n grows, and is eventually achieved when the total resummation is performed. It is interesting to notice 
that the "intercept" of g^ s DGLAP in this case could depend on Q 2 through the Q 2 -dependence of a^GLAP- Such a 
dependence originates from the Q 2 -dependence of a s . However, we have shown in Sect. IV that the parametrization 
a s = a s (Q 2 ) should not be used at small x. 
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